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Abstract 

We propose new textures for the fermion Yukawa matrices which are generalizations of 
the so-called Stech ansatz. We discuss how these textures can be realized in supersymmetric 
grand unified models with horizontal symmetry SU (3)h among the fermion generations. In 
this framework the mass and mixing hierarchy of fermions (including neutrinos) can emerge 
in a natural way. We emphasize the central role played by the SU (3)h adjoint Higgs field 
which reduces SU(3) H to U(2) H at the GUT scale. A complete 50(10) x SU(3) H model 
is presented in which the desired Yukawa textures can be obtained by symmetry reasons. 
The phenomenological implications of these textures are thoroughly investigated. Among 
various realistic possibilities for the Clebsch factors between the quark and lepton entries, 
we find three different solutions which provide excellent fits of the quark masses and CKM 
mixing angles. Interestingly, all these solutions predict the correct amount of CP violation 
via the CKM mechanism, and, in addition, lead to an appealing pattern of the neutrino 
masses and mixing angles. In particular, they all predict nearly maximal 23 mixing and 
small 12 mixing in the lepton sector, respectively in the range needed for the explanation 
of the atmospheric and solar neutrino anomaly. 
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1 Introduction 



In the last years the interest for the issue of neutrino mass generation has been renewed, 
especially by the SuperKamiokande data on atmospheric neutrinos These data can be 
explained by v„ — v T oscillation within the following parameter range (at 99 % CL): 

Sml tm = (1 - 9) x 10- 3 eV 2 , 

sin 2 26 atm > 0.8 (1) 

(the best-fit values are 5m 2 tm = 3 x 10~ 3 eV 2 and sin 2 2# atm = 0.99), while the explanation 
through the v^ — Ve oscillation is strongly disfavoured [§]. On the other hand, the solar neutrino 
data |J can be interpreted in terms of v e oscillation into u^, u T , or their mixture within the 
following parameter range j|, Bj: 

5m 2 ol = (3 - 10) x 10~ 6 eV 2 , 

sin 2 26» so i = (0.1 - 1.4) • 10" 2 (2) 

(the best-fit values are <5m 2 ol = 5 x 10 -6 eV 2 and sin 2 26* so i = 7 x 10~ 3 ), the so called small- 
mixing angle MSW solution ||.[] Therefore, the experimental data provide a strong evidence 
that like charged fermions, also neutrinos have masses and weak mixing. Namely, barring the 
less natural possibility that the neutrino mass eigenstates ^1,2,3 are strongly degenerate and 
assuming that 7713 > 7712 > m\, the 8m 2 ranges in eqs. ([l]) and (Q) directly translate into 

m 3 - (^atm) 1/2 = (3 - 10) • 10~ 2 eV, 

m 2 ~ (5m 2 ol ) 1/2 = (1.7 - 3.3) • 10~ 3 eV, (3) 

which indicates that the neutrino mass ratio 7713/7712 is similar to that of charged fermions, 
mr/rn^i or mb/m s . However, the neutrino mixing pattern inferred from (|l]) and @ is drasti- 
cally different from the well-established structure of the Cabibbo-Kobayashi-Maskawa (CKM) 
mixing of quarks. 

The explanation of the fermion flavour structure is one of the most challenging problems of 
particle physics, and the neutrino case constitutes a part of this issue. In the standard model 
the charged fermion masses emerge from the Yukawa terms :Q 

uiY^ qj H u + d^iqjHi + .;¥;'/,//,/. (4) 

where qi = (u,d)i, uf, d\ and b L = {v,e)i, e\ are the quark and lepton fields of three families 
{% = 1, 2, 3), and H u d are the Higgs doublets, with the vacuum expectation values (VEVs) v Ui d 
determining the electroweak scale, (tj 2 + v 2 ,) 1 / 2 = v\y = 174 GeV. There is no renormalizable 
term that gives rise to the neutrino masses. However, the latter can get Majorana masses from 
the lepton-number violating higher order operator 

li~Y l JljH 2 , Y u = Yj (5) 



3 We concentrate mainly on this solution, whose parameter range naturally emerges in the models presented 
below. One has to remark, however, that the solar neutrino data are still ambiguous and cannot discriminate 
other possible solutions. In particular, the recent analysis show that the large-mixing angle MSW and long 
wavelength "just-so" oscillation solutions have comparable statistical significance || p). 

4 In the following we mainly deal with the supersymmetric model, though the general features of our discussion 
are equally valid for the non-supersymmetric case. 
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which is cutoff by some large scale Ml, e.g. the grand unification or Planck scale (?]].[] Hence, 



the charged fermion masses are ~ vyy and the neutrino masses are ~ vfy/Mi which makes it 
easy to understand why the latter are so light. 

The quark mixing is originated from the misalignment between the Yukawa matrices Y u 
and Analogously, the neutrino mixing is due to misalignment between Y e and Y u . In 
particular, one can choose a basis where Y u and are diagonal, Y u = Diag(Y u , Y c , Y t ) and 
Y u = Diag(Yi, Y 2 , Y 3 ). In this basis, the quark and lepton mixing angles are determined 
respectively by the form of Y^ and Y e which in general remain non-diagonal. The latter are 
to be diagonalized by the bi-unitary transformations: 



U' d T Y d U d 

u: t ym p . 



rl) 



Diag(Y d , Y s , Y b ) 
Diag(Y e , Y M! Y T ), 



(6) 



with the unitary matrices U d , e rotating the left handed (LH) states and with U d rotating the 
right handed (RH) ones. Thus, in this basis the CKM matrix of weak mixing between the 
physical quark states (u, c, t) and (d, s, b) is V q = U d , while the matrix Vi = U\ determines the 
properties of the weak current between the charged leptons (e, /x, r) and the neutrinos. Namely, 
Vi relates the neutrino flavour eigenstates (f e , f M , u T ) to the mass eigenstates (1/1,1^2,^3), and 
describes the neutrino oscillation phenomena. 



The unitary matrices 



V„ 



^ V ud V us V u b 
Vcd Vcs Vcb 

\ v td V ts V t b 



( Vel V e2 V e3 



v, 



V, 



Ml 



Vu2 



(7) 



V v Tl v t2 v t3 



are usually parameterised as 



<i-i 



C12C13 



-S12C23 - Cl2S23«13e 

V S12S23 - c 12 c 23 s 13 e t5 



S12C13 si 3 e 
c\ 2 c 23 - si 2 s 23 s r3 e lS S23C13 
-C12S23 - si 2 c 23 si 3 e lS C23C13 



(8) 



g,i 



where Sy (cij) stand for the sines (cosines) of the three mixing angles 612, @ 23 and #13 while 
5 is the CP violating phase. (In the case of Majorana neutrinos, the lepton mixing matrix 
contains two additional phases 5\ y2 , factorised out as V\ ■ Diag(l, e iSl , e lS2 ), which, however, are 
not relevant for neutrino oscillations.) In the following, we distinguish the angles and phases 
in V q and VJ by the superscripts l q' and T, respectively. 

The explanation of the fermion mass and mixing pattern is beyond the capacities of the 
standard model. The Yukawa matrices Y Ujdje , are arbitrary complex 3x3 matrices, not 
constrained by any symmetry property. Concerning the neutrinos, apart from the Yukawa 
factors Y u in (|5|), also the lepton- number violation scale Ml is an arbitrary parameter. The 
magnitude of the latter can be inferred from the mass value m 3 in (||). Namely, from m 3 ~ 
Y 3 Vyy/ML we obtain 

Y^Ml ~ 10 15 GeV, (9) 

which is rather close to the GUT scale Mq ~ 10 16 GeV. In particular, the normalization of 
the operator (JsJ) as Ml = Mq implies Y 3 ~ 10. Let us remark that there is no contradiction 

5 Any known mechanism for the neutrino masses reduces to such an effective operator. E.g., in the 'see-saw' 
scheme M it is obtained after integrating out heavy-neutral fermions with Majorana masses ~ Ml- 
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2«d l^st 


yd j^nd 


Y c /Y u = 250 - 1000 


Y t /Y c = 230 - 330 (At=MwO , 320 - 1000 (M=Mf?) 


Y s /Y d = 17-25 


F&/Ys = 30 - 85^=^), 33 - 120 (m=Mg ) 


Y^/Y e = 207 


YtIYv = 17 


Y 2 /Y 1 > 1 


y 3 /y 2 = io-60 



Table 1: Quark and lepton mass ratios between second/first and third/second families. All mass ratios 
are independent of the renormalization scale p, except Y t /Y c and Y^/Yg (in the case of reasonably small 
tan/3 = v u /vd)- For definiteness, the latter ratios are evaluated at the electroweak scale, \x — Mw, 
while at the GUT scale, /i = Mq — 10 16 GeV, they are scaled respectively as -B t -3 and -B t _1 , with the 
factor B t — B(Y t ) varying from 0.9 to 0.7 for the top Yukawa constant at the GUT scale in the interval 
Y t G = 0.5 — 3 (see details in section 4). 



if the effective coupling constant Y3 is large. E.g. in the context of the see-saw mechanism |8|] 
I3 is determined by the ratio Y r ^ irac /Yp v iajorana which can easily happen to be ~ 10 even if the 
Yukawa constants loirac and YMajorana are order 1. 

One aspect of the flavour problem is related to the inter-family mass/Yukawa hierarchy 
among fermions |]lTj|]. For example, the top Yukawa constant is Yt ~ 1, while the first-generation 
constants Y e ^ u d are much smaller, ~ 10~ 5 . Concerning the neutrinos, they also seem to indicate 
a mass hierarchy similar to that of the charged fermions. In particular, eqs. (||) show that 
the ratio Y3/Y2 lies in the range 10 — 60. There are no data restricting the ratio Y2/Y1, but it 
may also be large. For the sake of comparison, the mass ratios between different families are 
summarized in the Table |l|. 

On the other hand, the mixing structure is very different between quarks and leptons (see 
Table §). The most striking feature is that the atmospheric neutrino data favour the nearly 
maximal 23 mixing of neutrinos, while the 23 mixing of quarks is very small. On the contrary, 
the MSW solution implies a very small 12 mixing angle for neutrinos as compared to the 
sizeable 12 angle of the Cabibbo mixing for quarks. 

It is widely believed that the fermion flavour structure can be understood within the con- 
text of the grand unified theories with horizontal symmetries. The latter acts between the 
fermion families and thus ca;ii help to constrain the form of the Yukawa matrices Y w ^ e v . In 
particular, it can dictate the hierarchical and closely aligned structures of Y u and and 
thus simultaneously explain the large mass splitting among different families as well as the 
smallness of the quark mixing angles. However, the origin of the large neutrino mixing is not 
a priori clear. Thus, the crucial point that must be explained by any realistic flavour theory 
is the dramatic difference between the quark and lepton mixing angles.^] 

In a previous paper [13] we have outlined that such a complementary pattern of quark and 
lepton mixings can naturally emerge in the context of grand unified theories. In the SU(5) 
model the quark and lepton states of each family are combined in the multiplets 5j = (d c , 
10j = (u c ,e c ,q)i, and the Higgs doublets H u and fit respectively into the representations 

6 In last years many efforts have been applied to derive the neutrino mixing pattern from particular textures 
of the neutrino mass matrix Differently from some point of view adopted in J!]]], we believe that a 

satisfactory scheme can be obtained in the context of a well-defined and complete flavour theory including the 
charged fermions. Such attempts have been taken e.g. in refs. |l2j where the neutrino mixing pattern was 
discussed in the context of grand unified and horizontal symmetries. 
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q(uark) 


l( epton) 


#12 


12.7 ±0.1 


2.4 ± 1.0 


#23 


2.3 ±0.1 


45 ± 13 


#13 


0.21 ±0.03 


< 15 



Table 2: Quark and leptonic mixing angles (in °). The values of 9 l 23 and 6[ 2 correspond to the 99% CL 
regions of the AN and SN data fitting, eqs. (0) and (|), while the upper bound for 0[ 3 comes from the 
combination of CHOOZ |l4| and atmospheric neutrino data [Q. The values of B\ 2 , &23 an d ^13 reflect 
the present data on |V^ S |, |V c fc| and | V^f, / V^t, | [^|. The last two values depend on the renormalization 
scale and for /i = Mq they should be scaled by a factor B t , i.e. become smaller by (10 — 30)% . 



H ~ 5 and H ~ 5. The terms responsible for the fermion masses read as: 

lOiG^'10^ + BiG^'lO.-ff + -l-EiG'jEjH 2 , (10) 

Ml 

where the Yukawa constant matrices G u and G v are symmetric due to SU(5) symmetry reasons 
while the form of G is not constrained. After the SU(5) symmetry breaking these terms reduce 
to the standard model terms (Q) and (||), with 

v C* v ( ' 

Y d = G, Y e = G T . (11) 

Hence, in the basis where the matrices G u and G u are diagonal, the rotation angles of the 
LH leptons coincide with that of the RH quarks: U e = U' d , and vice versa, = U' e . (One 
has to remark, that in grand unification the rotation angles of RH states also have a physical 
sense as they define the fermion mixing in B- and L-violating currents.) Hence, if G 23 ~ G 33 
while G 32 <C G 33 , the large 23 neutrino mixing will occur together with the small 23 mixing 



of quarks [12, 13]. In other words, the phenomenology requires that the matrix G should have 



a strongly asymmetric form. 

In particular, a very interesting situation emerges in the case of Fritzsch-like texture [15] 
with vanishing G 22 and G 11 entries. In addition, the entry G 33 can be taken as SU(5) singlet 
thus maintaining the b — r Yukawa unification, while the off-diagonal entries G 23 etc. should 
contain some SU(5) Clebsches to avoid the wrong relations Y^ s = Y e ^ of the minimal SU(5). 
This means that the masses of the third generation emerge from the Yukawa couplings to the 
Higgs 5-plet, while the lighter generation masses are contributed also by 45-plet. In fact, it is 
not necessary to introduce the elementary Higgs 45-plet. Instead, the off-diagonal entries G 13 
can be regarded as operators dependent on the adjoint Higgs $ ~ 24 of SU(5), G iJ = G 4 - 5 (<!>).[] 
In this case the relation Y^ = YjT is not exact anymore, but it fulfills with the precision of 
these Clebsch factors. Nevertheless, this property ensures that between the 23 mixing angles 



This should be understood as expansion series G y ($) = G (0)lJ + G (1)l3 |j + where M is some cutoff 
mass larger than the GUT scale (e.g. the string scale) . This means that the off-diagonal terms emerge from the 
effective higher-order operators j^HSiG^^Wj etc., just on the same footing as the last term in (|ici|). Since 
in general the operator $ • H is represented by the tensor product 24 x 5 = 5 + 45, it can distinguish the 
corresponding entries in the matrices Y e and Yd and hence avoid the lepton-quark degeneracy of the minimal 
SU{5) theory. 



4 



of quarks and leptons the following relation is fulfilled with a good accuracy [13]: 



tan #23 tan 



2.3 



m T mb 



1/2 



(12) 



This product rule indeed works remarkably well. It demonstrates a 'see-saw' correspondence 
between the lepton and quark mixing angles and tells us that whenever the neutrino mixing is 

23 



large, tan #23 ~ 1, the quark mixing angle comes out small and in the correct range, tan #23 



0.04. 



Namely, in ref. |T3[ 


we 




( Yu 





Y u = 


is 


Y c 
















Y u = 





Y 2 




V 






3| we have considered the following Yukawa pattern: 










' 


A d 


° ^ 


( 





A d 





H 


A' d 





B d = 




\A d 





B d 




V 


B' d 


D J 


\ 





\B d 


D 




/ 


A e 


^ 







A e 





Y T = 

1 e 







B e = 


( 


l A e 





B e 




V 


K 


D J 







I B e 


D 



(13) 



where the parameters b = B e /B' e = B d /B' d and a = A e /A' e = A d /A' d reflect the possible 
asymmetry between the off-diagonal entries in Y^ e . In addition, the Clebsches A e /A d = 
and B e /B d = ks are allowed to be different from 1 and their values can be extracted from the 
quark and lepton masses. 



Indeed, the texture predicts the following relations between the Yukawa eigenvalues |13| : 

Y s - Y d 1 



Y T 



Y„, - Ye 



k\z- 



Y d Y s 
YpY,, 



1 



k\z- 



and 



where 



m d 



2* -2 



k 2 

K A 
Zk A B 



m,, 



m, 



+ 



m, 



m. 



+ - k 



m. 



rrie 



B 



m 1 



(14) 



(15) 



(16) 



In particular, the phenomenologically correct picture for the fermion masses emerges if the 
asymmetry parameter b is large enough, b ~ 10, and when kj± ~ 1 and ks — 2 [|l^]. In this 



case the relation (15) leads to the correct prediction for the strange/down quark mass ratio. 
On the other hand, for the 23 mixing ang les in (@) we obtain the relations: 



tan #23 ~ b 



,V2 



tan#| 3 ~ b~ lf \ 



(17) 



from which the product rule (|l~2|) is immediate. These relations also point to a large asymmetry 
factor 6.0 In particular, for b ~ 10 the neutrino mixing becomes nearly maximal, # 2 3 ~ 45°, 
whereas the quark mixing angle gets small, #23 ~ (2 — 3)°. 

8 The symmetric case 6 = 1, i.e. Yd, e having the familiar Fritzsch texture, is obviously excluded, since the 23 
mixing of quarks is too big, #| 3 = (6 — 10)°, while it is too small for leptons, #23 = 13.7° (cfr. the experimental 
values in Table ^ . 



5 



In the general case, for arbitrary a, the product rule similar to that in ([l?]) can be obtained 
also for 12 mixing angles. However, from phenomenological grounds one can deduce that no 
significant asymmetry should occur in the 12 sector of the matrices (|l3|). Namely, the successful 
relation for the Cabibbo angle: 

tan 09, ~ x f^, (18) 



indicates that Ad ~ A' d , so that a = 1 is appropriate. Interestingly, in this case we obtain for 



the leptonic 12 mixing [13] 



1 I2 



/cos 9^. (19) 



Given that # 2 3 ~ 45°, this implies 9[ 2 ~ 3°, which is in the range relevant to the MSW solution 
of the SN problem (cfr. Table |). 

Therefore, for a = 1 the ansatz ( |l3| ) depends on six parameters, so it is highly predictive. 
Among these, three Yukawa entries A e ,B e ,D can be expressed in terms of lepton couplings 
Ye,fi,r and b. Thus 9 physical quantities, the Yukawa eigenvalues Yd )S ^ and the mixing angles 
in V q and Vi, are left as functions of three Clebsch factors, ks,kA and b. Hence, at the GUT 
scale these are connected by six relations like ([il]), (|l7|), etc. As discussed above, in the first 
approximation these relations well reproduce the observed pattern of the quark and lepton 
mixing angles for b ~ 10 and ks — 2, kj± — 1- However, a closer inspection shows that for 
precise Clebsch values ks = 2 and k^ = 1, which could be motivated in the GUT context, 
and a = 1 which could follow from the U{2)h horizontal symmetry, the predictions for \V US \ 
and I V c b\ substantially deviate from the corresponding experimental values which presently are 
known with a very good accuracy. The case summarized above implies a real CKM matrix 
since Y u is taken diagonal and so all phases in Y^ can be absorbed by the field re-definitions. 
Nevertheless, it would be interesting to address the CP issue in a scenario in which Y u is still 
kept diagonal. 

In the present paper we suggest a new grand unified ansatz for the fermion masses, which 
is an extension of the pattern (|l3|) of ref. [13|. Namely, we consider the case when the 



matrices G u and G v are both diagonal and related as G u oc b /_1 • G u , with b' being a 
diagonal matrix b' = Diag(l, 1, b'). Regarding G, we assume that it has the diagonal part 
oc G u and the off-diagonal one oc b" 1 • A($), where A is an antisymmetric matrix with <3?-field 
dependent entries inducing different Clebsch factors for the down quarks and charged leptons, 
and b = Diag(l, 1, b). In other words, we consider the Yukawa textures: 

Y u = Yg, Y d = P Y U + b^Aa, 
Y v = r/b'^Yu, Yj = P Y U + b- 1 A e , (20) 

where p and rj are proportionality coefficients. For 6 = 1 this pattern resembles the so-called 



Stech texture proposed long time ago by Stech |16|] and independently by Chkareuli and one 



of the authors [17], where the off-diagonal entries of Y^ were assumed to be antisymmetric, 
Yd = pY u + Ad- This case, however, is completely excluded by experimental data. As for 
our texture, the matrices b _1 A^ )e remain antisymmetric only in the 12 block, while the other 
entries are skew due to the factor b 7^ 1. In the explicit form the Yukawa matrices read as: 

/ Yu 
Y u = Y c 
\ 




6 



Y„ 



( Yu \ I \B 

o y c o I , y t 



o o 





A e £D B e I . (21) 



Clearly this texture represents an extension of the Fritzsch-like 'zero-texture' (|13| ) considered 
in [13]. In particular, for vanishing 13 and 31 entries, Cd, e = 0, this texture essentially reduces 



to the latter as the diagonal 11 and 22 entries in ^Yd,e ( |20D are quite small, since the ratios 
Y u /Yt and Y c /Yt are much smaller than 1/3 and 2/3 Yukawa ratios in the down quark and 
charged leptons. In practice, as we show below, one can safely take the limit Y u /Yt — > and 
ignore the 11 entries in Y^ e . As for the 22 entry, though it is small, it can provide significant 
corrections (in particular, to V^). In principle, it also contains the relative (un-removable) CP 
violating phase which however cannot provide a sufficient amount of CP violation. This can 
be remedied by non- vanishing C e ^- 

Moreover, we show that for large values of b (~ 10), the texture (|2l| ) provides a very 
appealing description of the quark and lepton mass and mixing structures. Interestingly, 
the case b' ~ b ~ 10 also provides nice relations between neutrino and up-quark masses, 
Y 2 /Y 3 = b%/Y t (cfr. Table 1). 



Hence, the complete pattern (21) is first motivated on the phenomenological point of view 
as it can fit the physical observables with excellent precision. From a theoretical point of view, 
we show that these textures can be justified in the framework of grand unification by invoking 



also the concept of horizontal symmetry among three fermion families |L^, 17]. 

The paper is organized as follows. In Sect. 2 we discuss some theoretical tools for building- 
up the textures (p0|) in a realistic and predictive manner in the framework of supersymmetric 
SU(5) model with SU(3)h horizontal symmetry. The discussion in Sect. 3 is meant to provide 
an existence proof of those realizations in the framework of S'O(IO) x SU(3)h model. Three 
particular ansatze are singled out, characterized by different assignments of the Clebsch fac- 
tors between the quark and lepton entries, and in Sect. 4 their phenomenological analysis is 
presented and discussed in all details. Finally, we summarise our findings and conclude. 



2 Theoretical framework: GUT and horizontal symmetry 

In this section we would like to sketch the general ideas needed to obtain the Yukawa textures 



( P0[ ) in the SU(5) model with the SU(3)h horizontal symmetry [18, 17|. In the course of the 
presentation, it will become apparent that the discussion has to take a less general character 
as for example some additional symmetries should be invoked. This aspect will be faced in 
Sect. 2.2. In the discussion SU(5) is taken as a prototype theory featuring the basic properties 
of grand unification. In a more general context one could think of GUTs based on larger 
groups (e.g. SU(N) with N > 6) which contain SU(5) and perhaps also unify both SU(5) and 
SU(3)h- A particular model based on S'O(IO) x SU (3)h will be presented in the next section. 

2.1 General aspects 

We consider the supersymmetric SU(5) model with the horizontal symmetry SU(3)h where 
three fermion families are unified in chiral superfields: 

5; = (d c , l)i ~ (5, 3), 10i = (u c , q, e c )i ~ (10, 3), (22) 
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(i = 1,2,3 is SU(3)h index), while the Higgs superfields 5,5 are singlets of SU(3)h- 

F~(5,l), ff~(5,l). (23) 

We also assume that the theory is invariant under R-parity. In other words, we impose the 
matter parity Z 2 under which the matter superfields change the sign while the Higgs ones are 
invariant.^ 

Since the fermion bilinears transform as 3x3 = 3 + 6, their "standard" Yukawa couplings 
to the Higgses are forbidden by the horizontal symmetry.^ Hence, the fermion masses can be 
induced only by higher order operators involving a set of "horizontal" Higgs superfields X 13 
in two-index representations of SU(3)h- symmetric X l J = S**^ ~ (1,6) and antisymmetric 
X$ = Afel = £ ijk A k ~ (1,3):[] 



M M J M 2 3 



where M is some large scale (the flavour scale). Needless to say, by SU(5) symmetry reasons, 
the antisymmetric Higgses A can participate only in O. Therefore, in the low-energy limit the 
operators (|2~4]) reduce to the Yukawa couplings: 

Y "~¥ ' Yd ' Ye m ' W L ~W (25) 

In this way, the fermion mass hierarchy can be naturally linked to the hierarchy of the 



horizontal-symmetry breaking scales [17, 20 1. For more details on the horizontal VEV struc- 
tures, see e.g. pl| . 

In particular, let us assume that the horizontal Higgses include one sextet S with a VEV 
taken diagonal: 

/ Si \ 

(S ij )= S 2 , 5 3 »5 2 »5i ! (26) 
\ S 3 J 

and a set of triplets A^ = e ljk A n j t , n = 1,2,3, which in general have the VEVs towards all 
three components: 



V 




Ax > A 2 > A 3 - (27) 



(see Appendix A). Then from (25) we see that Yf ~ 1 implies £3 ~ M which indeed can 



naturally arise from the Higgs sector as we show in Appendix A. Similarly one can expect 

9 The SU(3)h symmetry itself does not not prevent the i?-parity violating couplings e^^WSjEk from where 
the B — L violating terms u c d c d c ,qd c l, e c ll arise. For a general discussion on the relation between R parity 
and horizontal symmetries, see [n9[. In particular, in case of the horizontal symmetry SU(4)h the R-parity 
could emerge as an accidental symmetry. One has to remark, however, that in many cases considered below the 
matter parity automatically follows from the discrete symmetries imposed on the model. 

10 The following discussion can equally apply in both the cases of global and local horizontal symmetry, 
though it would be more appealing to regard SU(3)h as gauge symmetry, like 5(7(5). 

11 The theory may also contain conjugated Higgses Xij in representations S ~ (1,6) and A ~ (1,3). These 
usually are needed for writing non-trivial superpotential terms in order to generate the horizontal VEVs (see 
Appendix A). The conjugated fields, however, do not couple the fermion superfields (^2|) and thus do not 
contribute to the quark and lepton masses. 
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that also Y\, T ~ 1 which would require large tan (3 regime. However, as we shall see below, in 
realistic schemes also small tan/3 can be naturally accommodated. In addition, the operator 
O u translates into Y^/Ml ~ S3/M 2 ~ Yt/M. Therefore, the flavour scale M is of the order of 
the B-L violating scale Ml ~ 10 15 GeV (cfr. (^)). The closeness of the flavour scale with the 
GUT scale Mq ~ 10 16 GeV suggests that they may have a common origin. (The mismatch 
between the estimate of M and Mq can be due to some spread of the coupling constants in 
the theory.) 

It is well known that the effective operators (^) may arise entirely from a renormalizable 
theory, after integrating out some heavy degrees of freedom |2^] , and thereby the flavour scale 
can be related to the mass scale of the latter. It is also plausible to think that this mass scale 
is set by large (order GUT scale) VEV of some Higgs O, ~ M. 

In particular, one can introduce vector-like states in the 10- and 10-representations of 
SU{5): 

T = (U c , Q, E c y ~ (10, 3) , T t = (U, Q c , E) t ~ (ID, 3), (28) 
and the SU (5)-singlet fields 

iV*~(l,3), ^-(1,3), (29) 
and consider the following Yukawa couplings in the superpotential: 

W T = WHT + 5HT + TX10 + TttT 

W N = 5HN + X S N 2 + NtiN (30) 

(family indices are suppressed, and order 1 coupling constants are understood in each term). 
Clearly, in Wt the horizontal Higgs X can be both of the type X s = S or X a = A, whereas 
in Wjsr only the symmetric Higgses X s can participate. Since T, T contain heavy states with 
electric charges 2/3, —1/3 and —1, their exchange can induce the operators O u and O relevant 
for the masses of all charged fermions: up quarks, down quarks and leptons. Analogously, the 
operator O v for the neutrino masses emerges via the exchange of the singlet states N, N. The 
relevant diagrams are shown in Fig. |l[ 

Therefore, after integrating out the heavy states, we obtain the operators ( pl|) in the form 

O u = WiHM^X + X t Mt 1 H)W = 2 • lOff (M^ 1 5)10 + lOi^M^ - AM^)W, 
O = 5(FMy 1 X)10 = BHiM^S + M^ 1 yl)10, 

o v = mM^SM^m (31) 

where My, Mat are 3x3 mass matrices of the heavy states, induced from their couplings to 
the Higgs fi. Therefore, if £1 is a gauge singlet, then My, Mat will be degenerate in family 
space, M/r,jV ~ M • 1, where 1 = Diag(l, 1, 1) is the unit matrix of SU(3)h- In this case, if X 
are SU(5) singlets, only the symmetric fields X s = Xj can contribute O u in ([H]), while the 
contributions of the antisymmetric ones X a = — X^ cancel out. On the other hand, both X s 
and X a can contribute to O. Therefore, the operators ( |3l|) reduce to the Yukawa couplings of 
the form: 

Y u ~ My 1 (5), Y d ,Y? ~W T l (S)+W T l (A), Mj}Y v ~ M Af 2 (5). (32) 
More explicitly, these can be written as: 

Y u tx S, Y„ oc S, Y d = Yj oc S + A, (33) 
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Figure 1: The heavy fermion exchanges giving rise to the effective higher order operators for the 
fermion masses: (a) for up quarks, (b) for down quarks and charged leptons, (c) for neutrinos. 



where S = M~ 1 (S) and A = M~ 1 (A) are respectively the symmetric and antisymmetric 
matrices reflecting the form of the VEVs ( |26| ) and (p7|). (We omit the dimensionless factors 
which keep track of the coupling constants in (|30|) and are not specified at the moment). 

It should be immediately noted, however, that the texture ( |33| ) is far from being realistic 
for many reasons: 

(A) Yd oc Y u + A with antisymmetric A leads to very big 23 mixing, — (Y^/Yf,) 1 / 2 ~ 
0.1 — 0.2, much above the experimental value of 

(B) the relation ~Y U oc Y u implies Y\ : Yi : I3 = Y u : Y c : Yt and thus too big a hierarchy 
between the neutrino mass eigenvalues and 7712 (cfr. Table HJ). 

(C) finally, Y^ = Yj features the minimal SU(5) degeneracy of the Yukawa couplings, Y^ s ,b = 
Y e „ T , whose drawbacks have been already recalled. 

The latter problem can be easily cured by making use of the adjoint superfield <I> which 
breaks SU(5) down to SU(3) x SU(2) x U(l) by the VEV proportional to the hypercharge 
generator: ($) oc \ y = Diag(l/3, 1/3, 1/3, -1/2, -1/2). Thus if the matrix A in © effectively 
contains $>, A = A(<&), effective Clebsch factors emerge discriminating quarks from leptons. 
Therefore, at least some of the antisymmetric Higgses A can be taken in the mixed representa- 
tion A ~ (24, 3).0 The tensor product is 5 x 24 = 5 + 45 + 70, where only the first two terms 
are relevant for the fermion bilinear 5 x 10 = 5 + 45. Then the unwanted relations Yd s = Y e ^ 
can be avoided while maintaining the successful relation Yj, ~ Y T . 

The first two problems can be solved in a similar way, by making use of the adjoint Higgs 
of SU(3)h- The masses of the heavy fermions transform as 3x3 = 1 + 8 representations. 
Hence one can assume that the set of Higgses £1 contains, besides the singlet / ~ (1,1), also 

12 Alternatively, one can take all the antisymmetric Higgses as 577(5) singlets, A ~ (1,3), but assume that 
some of them act in combination with <E> ~ (24, 1), to end up with the effective operator |j|5 • 10H (cfr. (p4|)) 
so that the tensor product H ■ $ ■ A effectively contains both (5, 3) and (45, 3) representations. Such an attitude 
will be taken in constructing the 5*0(10) model in Sect. 3. 
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the superfield £ ~ (1, 8) with the VEV (£) pointing towards the As generator :P[ 




(£) ~ M ■ 1 . (34) 



This VEV breaks SU(3)h down to SU(2)h x U(1)h and thus discriminates the third fermion 
family from the others ^3|. Hence, the heavy-fermion mass matrices Mjyv maintain only the 
SU(2)h invariant form M = M ■ Diag(l, 1, b), with b in general different from 1. In this way, 
from (^) we can obtain the desired Yukawa pattern (pC|).P| 

Notice however, that for achieving our goal, we have implemented the adjoint structures 
(S into M and <3? into A) in eq. (^2|) which is not the right place. Indeed, (|32| ) follows from 
( |3"l"D if and only if neither M nor A contain adjoint representations of SU(5) or SU(3)h- For 
example, if A contains the SU(5) adjoint, then the tensor product A x H ~ (24,3) x (5,1) 
effectively contains the (45,3) representation which would induce off-diagonal antisymmetric 
contribution also to Y u . On the other hand, if ^4's are SU(5) singlet, but My contains the 
SU(3)h octet, then (1,3) x (1,8) x (5,1) effectively contains (5,6) which would still give rise 
to symmetric off-diagonal entries in Y w . Therefore, in either case the antisymmetric VEVs 
(A) contribute the first term in (|24|) and so induce off-diagonal entries in Y u ruining in such 
a way the desired pattern. Tackling this issue requires more theoretical ingredients and tools 
which will be discussed below. 



2.2 Some specific realizations 

Thus, we have seen that in order to arrive at the desired textures (|20|), both the SU(3)h 
singlet / ~ (1, 1) and octet £ ~ (1,8) are needed for generating the heavy fermion masses. In 
addition, the diagrams contributing via the S and A type Higgses are to be "differentiated" 
by some additional symmetry so that S will contribute to all the matrices Y u , Y^ and Y^ e , 
while A's only to Yd, e - Hence, we have to assume that besides the local SU (5) x SU (3)h sym- 
metry, the theory is invariant under some set Q of abelian symmetries which may contain e.g. 
non-anomalous or pseudo-anomalous gauge U(l) factors, continuous or discrete R-symmetry, 
discrete symmetries like Z2, Z3, etc. 

First, on the basis of the considerations outlined in the previous section, we put forward the 
needed couplings in the superpotential. Second, those will be motivated by some additional 
symmetry of Q. 

Now about the first step. We assume that the flavour scale is determined by the VEVs of the 
singlet / ~ (1, 1) and the octet S ~ (1,8), respectively (I) = M and (£) = xM-Diag(l, 1, -2), 
x ~ 1. Then the combination I + X provides order M masses to the heavy fermions, splitting 
the third generation from the first two and maintaining the latter degenerate. As in the 
previous subsection, we take the horizontal VEVs in the form ( p6| ) and (^7|) and define the 
dimensionless VEV matrices as S = M^ 1 (S) and A n = M~ 1 (A n ). 

13 In the following of the paper, the symbol f2 stands for a set of Higgs fields in real representations of the 
gauge groups. When needed the content of Q will be specified as singlet I ~ (1, 1), octet E ~ (1, 8) or 24-plet 
$ ~ (24, 1). 

14 In group-theoretical language this can be rephrased as follows. The tensor product A ■ E reads as 3 x 8 = 
3 + 6+15 where both 3 and 6 can match the fermion bilinears 3x3. In this way, the off-diagonal entries in ~Y e ,d 
are not anymore antisymmetric - only the 12 sector keeps on the antisymmetry owing to the residual SU(2)h 
symmetry. 
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• We further assume that the superpotential Wt (1^) related to heavy 10-plets (^8|) involves 
some set of horizontal Higgses X, which in fact consists of the sextet S ~ (1,6). Therefore, 
Wt becomes: 

W T = flOHT + glHT + X T TXW + T(a T I + #rE)T, (35) 

where we explicitly indicate the order 1 coupling constants. So the diagonal entries in Y M and 
Yd, e are induced respectively via the diagrams (a,b) of Fig. |l[ In order to respect the b — r 
Yukawa unification, the 24-plet $ ~ (24, 1) should not couple to TT. 

• The couplings in Wn are left as in ([3(]) and they induce ~Y V via the diagram (c) in Fig. 

i _ 

W N = h5HN + X N SN +N(a N I + (3 N E)N. (36) 



Thus, Y u is diagonal in the VEV basis (26). 

• Let us consider another set of horizontal Higgses X' , which contains antisymmetric fields 
A n in the representations (1,3) and (24,3), and assume that they contribute only to Y<j e via- 
the exchange of other heavy fermions. The simplest possibility is to introduce additional states 
in 5- and 5-representations of SU(5): 

T = (D c , Lf ~ (5, 3), Ft = (D, L c )i ~ (5, 3) , (37) 

having the following couplings in the superpotential: 

W F = X F 5X'F + g'FH'W + F{a F I + F E)F, (38) 

where H' is another Higgs 5-plet (or 45-plet) of SU(5). The corresponding diagram is shown 
in Fig. §(a). In the following, such a scenario will be called F-scheme. 

Alternatively, instead of the F + F one can involve an additional pair of 10-plets 

T H = (U' c , Q\ E'J ~ (10, 3) , % = (W, Q' c , E')i ~ (TO, 3), (39) 

with the superpotential terms: 

W T > = g'lH'T' + X' T T'X'W + T'{a' T I + (3' T T,)T'. (40) 

The corresponding exchange is shown in Fig. ^(b). This scenario will be referred to as the 
T'-scheme. 

In either scheme, the standard Higgs doublet should be regarded a superposition of the 
doublet components in H and H': = — s w H2, which is rendered light after arranging 

the doublet-triplet splitting, while the orthogonal combination s w H' 2 + c u H2 is left heavy, with 
mass ~ Mq (here a; is a mixing angle and c^s^) = cosa;(sina;)). On the other hand, we 
assume that H u comes entirely from H, while the field H', though present, is just a spectator 
in the Yukawa sector. On the contrary, Hd is contained in both H and H' with weights and 

respectively. 

Once the fields A n = (X a ) n contain the mixed representation (24,3) and/or H' is 45-plet, 
the diagrams in Fig. [2] induce off-diagonal entries in Y^ and Yj with different Clebsch factors. 
Let us consider in some more details e.g. the F-scheme. After substituting the VEVs of the 



horizontal Higgses, the couplings (|35|) , (|36|) and ( |38|) give rise to the following field-dependent 
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Figure 2: The diagrams inducing the off-diagonal contributions to ~Yd,e'- (a) exchange of 5-plets (F- 
scheme), (b) exchange of additional 10-plcts (T'-scheme). 



mass matrices for the up quark and neutrino states: 
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and for the down quark and charged lepton states: 
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(42) 



where T\^ e are the Clebsch factors dependent on the SU(5) content of H'. Namely, r/ e = 
if H' is 5-plet and r\ e = —3% if H' is 45-plet. In the above, each entry is a 3 x 3 matrix in 
the SU(3)h space. In particular f = /l, h = hi, etc. are flavour-blind (SU(3)h degenerate) 
matrices. 

As long as the heavy fermion mass matrices are induced by I and £, they will be diagonal 
and have SU(2)h invariant shape = o$Mb^ (\& = T,F,N), where 




a^p = acq, + xf3q>, 



(43) 



These will give rise to the deformation factors b and b' in fl20|) whose relations with the factors 
b^, are easy to catch - at any rate they will be shown later. 

Thus, the main information on the fermion flavour pattern is contained in the matrices 
^hj,,d,l,v The matrices m u ^ u are proportional to the VEV (5) = Diag(<Si, S%, S3), m u = Ay MS 
and nij, = X^MS. As for the matrices m^ e , they are antisymmetric and emerge from the 
VEVs of the triplets A n sketched in eq. (^7j). As far as the latter set contain (1,3) and (24, 3) 
representations, all those entries in general are non-trivial 5x5 matrices, Ak = -4.fc(U + Ck^y), 
where 11 and X y are the unit and hypercharge matrices in SU(5) space, respectively and £fc 
measures the relative weight of the singlet and 24-plet VEVs at each component (k = 1, 2, 3). 
Thus, the corresponding entries in and m^" differ by Clebsch factors C^' e = 1 + CkVd,u 
where = 1/3 and yi = —1/2 are the hypercharges of the d c and I fragments of 5-plets. 
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Therefore, after integrating out the heavy states, we obtain an effective low energy theory 
with the following Yukawa matrices: 



Y u = /(M r 1 m„ + m^M^ 1 ) 
Y d = gSuM^nLu + g' c^diiidMp 1 
= gSuM^niu + g' c^miMp 1 



and 



h 2 M N 2 m u . 



h 2 X 



N 



Therefore, Y u and Y„ get contribution only from S: 

a T T ' Ml a%M 

and thus both have diagonal forms related as: 

V Y 



where b' = Diag(l, 1, b') , and 



n 



h Xnclt 



(44) 
(45) 

(46) 
(47) 
(48) 



As for Y ej rf, they have a diagonal contribution gs u /(2f)Y u from the first term in ([44]), and 
off-diagonal entries A^b^ 1 from the the second term. Therefore, the full matrices Y^ e can be 
rewritten as 

Y d = P Y U + b- 1 A d , Y e r = pY u + b- 1 A e , (49) 

where p = ^js^ and b = Diag(l, 1, b), with b = bp , while the other dimensionless coefficients, 
such as the Clebsch factors related to the representations of H' and A, are absorbed into the 
definition of the entries of the antisymmetric matrices 



(50) 

Putting everything 
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so that the relative Clebsch factors are defined as ki 
together, we can write the Yukawa matrices as: 
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(51) 



The following comments are in order. 
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• First of all, Y% ~ 1 implies that S3 ~ brM, while the relations Y u : Y c : Y t ~ e 2 : e : fr^ 1 , 
well fits the upper quark mass hierarchy for e ~ 10~ 2 and bx ~ 0.1 — 1. As already noted 
in Sect. 2.1 this situation can naturally emerge from the pattern of the horizontal Higgs 
superpotential (see Appendix A). 

• Since the (3,3) element in Yd, e is also related to largest scale S3, we have approximate 
b — t Yukawa unification at the GUT scale: Y T ~ Y b ~ pYf. On the other hand, as in principle 
the mixing angle u between the doublets in H and H 1 can be small, i.e. sinu; < 1, the 
case Yb/Yt <C 1 is plausible even if the constants / and g are both of order 1. Therefore, the 
regime of moderate tan (5 is quite a natural possibility, in which case the renormalization group 
equations can be substantially simplified. In the following we assume that regime, though the 
case with large tan (5 should not be excluded. 

• The large 23 lepton mixing implies A\ ~ Y T . This also seems natural in the context of 
the model as the largest triplet VEV A\ should be of the order of the flavour scale scale M 
(see Appendix A). On the other hand, the small 23 mixing of quarks needs Af/b < O.lYf,, i.e. a 
large asymmetry parameter, b ~ 10, which in turn implies bp ~ 0.1. In this case we observe an 
inverse hierarchy between the first two and third family of the heavy F states, i.e. Mf 2 > M% . 

• Concerning the neutrino Yukawa eigenvalues, they obey the relations 

Y = Y <<h % = h Y (52) 

which also points to large b' (cfr. Table [l]), b' ~ 5 — 100. In particular, the relation b' = b 
could work, which can occur if bp = bp and b^ = 1 - that is when the same singlet-octet 
combination I + £ acts in Wt and Wf, while in Wm only the singlet / is present (/3at = in 
©)• 

• On the other hand, we obtain (cfr. 



ML = „ 10 i 5 GeV (53) 

Y3 Y t 

Since Y t ~ 1 and b' ~ 10, this translates into 7] ~ 10 2 if M ~ 10 16 GeV, i.e. order GUT scale. 
It should not be, however, very surprising to find 77 ~ 10 2 due to some spread of the coupling 
constants in the theory (cfr. 

Let us remark that the set X appearing in eq. ( |35| ) could contain also a triplet Higgs 
^3 ~ (1,3) having a VEV towards the third component, (A^) = e^^A^- This would induce 
antisymmetric 12 entries in the matrix m u . This, however, would not alter the diagonal form of 
Y n , since in the combination M^ 1 m M + m^M^ 1 the antisymmetric entries cancel out (thanks 
to the residual SU(2)h symmetry in My). However, in Y^ e this will induce antisymmetric 12 
entries with Clebsch factors k e = k^. 

The T'-scheme can be developed along the same lines. The only difference will be that in 
this case we have b = bj", so that there should be a direct hierarchy between the first two and 
third family of the heavy T 1 states, i.e. Mf 2 < M% . In this case b' = b can be obtained if 
br' = bjsr = bp, i.e. when Wtt' n share the same / + X combination. Alternatively, when 
6jv = 6v but bx = 1 (in Wt contributes only the singlet /), we would obtain another interesting 
possibility b' = b 2 which can be also compatible with the experimental data. Namely, it would 
point towards smaller 3/2 hierarchy between the neutrino mass eigenvalues: Y3/Y2 ~ 10 (cfr. 
Table 0). 
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We now come to the second step to prove that the superpotentials (j35|), (|36|) and (|38|) 
can be justified by some extra symmetry Q. First of all, we assume that Q contains a discrete 
symmetry 1Z under which all superfields in the theory as well as the superpotential change sign. 
This will allow only trilinear terms in the superpotential (provided that they are consistent 
also with SU(5) x SU(3)h and other symmetries in Q), while some singlet Y can have also 
linear term h?Y . In this situation all VEVs (including the GUT and horizontal ones) can 
emerge from a single mass scale A. 

Second, we introduce a Z 3 symmetry acting on the superfields as tp — > ip exp(^zQ^), with 
the following 'charges' Q^: 

Q = 0: Y, H'; T, 

Q = l: I,Z,S, A n , H, H'; 5, F, F, N, 

Q = -1 : I, S, An, H; 10, T, N. (54) 

where as agreed $ ~ (24, 1), E ~ (1,8) and /, / and Y are some singlets, ~ (1, 1), 

The most general superpotential invariant under 7Z x Z% is restricted to have the following 
form: 

Wkiggs = A 2 ^ + y 3 + Y$ 2 + ^ 3 + YII + I 3 + I 3 + IT? + S 3 + 
Y{SS + A n A n ) + S 3 + S 3 + A ± A 2 A 3 + A 1 A 2 A 3 + 

h{y + + i?'(y + §)ti + fr'/ja"' (55) 

and 

HVuk = 10(.ff + H')T + 5OT + r(J + S)T + T510 + 
+ F(J + S)F + FH'VQ + 
5FiV + NYN + SA^ 2 (56) 

respectively in the Higgs and Yukawa sectors. 

It is worth observing the following features of the Higgs potential. It has only one mass scale 
(in the linear term) which determines the physical scales present in the theory (GUT, flavour, 
etc.). Notice that among other degenerate minima, this superpotential allows a solution when 
all singlets, adjoints and sextet/triplets get VEVs linked to the scale A, modulo unspecified 
coupling constants in (^). For example, the superpotential terms for S and S considered in 
Appendix A, emerge from here with fj, ~ (Y) . 

There are no couplings between S, A and E, so the superpotential does not fix the rel- 
ative orientation of their VEVs. However, that will be fixed by supersymmetry breaking 
terms, along the lines discussed in the Appendix A. For example, these can be D-terms like 

/ d9 4 zz[aS^T, jf T,S + X n Al l T, ji 'EA n ] etc., where z = m6 2 stands for supersymmetry break- 
ing spurion. Depending on the form and sign of the coupling constants, these terms can fix 
the relative VEV orientation of S, S and A n in the SU (3)h space. In particular, the choice 
Ai 2 > favours "orthogonal" VEV directions of A\ 2 and S, while a < prefers the "parallel" 
VEV directions of S and E. Thus, in this case the largest entry in S is positioned at the same 
place as in E. The same for other fields. A typical consequence of the fact that the horizon- 
tal VEV configuration is dynamically fixed by small terms, is the presence of light horizontal 
Pseudo-Goldstone Higgses in the theory, reminiscent of the familons, even if the horizontal 
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symmetry is local As for the SU(5) part, the adjoint $ gets a VEV ~ A and, as far as 
the superpotential contains also <3? 3 term, no light fragments of $ are left behind. 



Notice that the term HIH' is not present in (55) although it seems to be allowed by the 
symmetry. The point is that if H' is a 5-plet, then as it has the same charges as H , that term 
can be rotated away by redefinition of H and H' . Instead, if H' and H' are 45-plets, then that 
term is forbidden by the SU(5) symmetry. Thus, after substituting the VEVs, the bilinears of 
the doublet and triplet fragments in H, H' and H , H' have the form: 



XiY + A 2 $ 

A 3 y + a 4 $ 



(57) 



where Ai, A2, • • • are 0(1) coupling constants previously understood. For rendering the Higgses 
H u d light, the 11 element of (|57| ) in the corresponding doublet sector has to be fine tuned by 
canceling big contributions between \{Y and A2 < 3?. (Then the triplet components will remain 
heavy.) Therefore, H u is contained entirely in H whereas emerges as a superposition of 

the doublet fragments in H and H', = —s u H2 + c^H' 2 , with tanw ~ (A 3 y+A 4 <s>) > wriUe the 
other combination of these doublets is superheavy. 

The Yukawa superpotential (|5^) contains all terms of eqs. ([35|), ( ptf ) and ( |38| ) which lead to 
the Yukawa pattern (|49|). Furthermore, some specific features have arisen. The heavy-states 
N mass matrix emerges only by the VEV of singlet Y, so bjy = 1 and thus b' = Kj} . Then, 
if I + £ acts in the same reducible combination in Wt and Wf, i.e. &t = then b' = b, 
which is phenomenologically welcome (52). For example, this could be the case if the reducible 



combination / + £ emerges from an adjoint Higgs of a larger group containing SU(3)h, e.g. 
SU(N + 3)h- Following similar lines of reasoning, also the T' scheme can be justified. 

Finally, in either schemes, the Clebsch coefficients k n = A^/A^, n = 1,2, 3, depend on the 
SU(5) content of the tensor products A n ■ H' . Needless to say that whatever values of the 
Clebsch factors can occur when the Higgses A n are in ~ (1 + 24, 3) with unspecified weights £ n 
of 1 and 24 component. However, it would be more attractive (though it as a matter of taste) 
to find a more "familiar" origin of the Clebsch factors phenomenologically needed. In other 



words, we assume (on similar grounds as e.g. in refs. |24]) that different entries in the Yukawa 
matrices are differentiated by some "clean" Clebsch factor k n in which the VEV structure 
of the original GUT is encoded through a specific representation of the intermediate heavy 
fermions (e.g. F or T' in our case). Therefore, let us think that at each entry, the Clebsches 
k n can have some specific rational values such as 0, ±|, ±|, ±1, ±|, ±2, ±3, etc. , possibly 
related to various charges (electric, hypercharge, B — L etc.). 

In particular, in the context of SU(5), whenever the product A n ■ H' is in 5- or 45-channel 
(e.g. simply when A is singlet, A ~ (1, 3), and the Higgs H' is in 5 or 45 representations), we 
have k n = 1 and k n = —3, respectively. In Table |3| we report the values for k n attainable for 
A ~ (1,3) and A ~ (24,3) (or M ~ 24), too. 

As in the next we shall consider in detail specific Clebsch patterns for the Yukawa matrices, 
characterized by different values of k n , now we would like to keep the discussion on the Clebsch 
factors more general to show how we arrive at those specific choices. 

Indeed, we could think of all these different ansatze in terms of larger vertical and/or 
horizontal symmetries. Namely, in these cases SU(5) may emerge just a subgroup of some 
larger GUT symmetry group. (The SO (10) example will be considered in the next section). 
Then the adjoint representation would generally emerge accompanied by a singlet partner. 
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SU (5)rep. 


A~ 1, 


,4 ~ 24, 


A~ 1, 






M ~ 1 


Af ~ 1 


M ~ 24 


F — scheme 


H' ~ 5 


1 


-3/2 


-2/3 




i7'~45 


-3 


9/2 


2 


T' — scheme 


H' ~ 5 


1 


6 


1/6 




i7'~45 


-3 


-18 


-1/2 



Table 3: Possible values for the Clebsch factors k n — A^/A^, as they can emerge in the F- and 
T'-versions for various assignments of the SU(5) representations to the Higgses H' and A and to the 
heavy-state masses M. The case when both A, M ~ 24 gives the same pattern as A, M ~ 1. Clearly, 
in each case the Clebsches obtained for 5' ~ 5 and ff' ~ 45 differ just by a factor —3. 



Hence, the corresponding combination 1 + 24 can have a VEV with two different eigenvalues 

($1+24) = Diag(x, x, x, y, y) (58) 



In the context of F-scheme (exchange of heavy 5-plets) this would produce a Clebsch k 
In the case of 10-plet exchange (T' scheme) the Clebsch value would be k = -prr- 



x+y ' 

For example, in the context of the SU[p) model, which, as a matter of fact, provides a natu- 



ral doublet-triplet splitting, the adjoint Higgs (35-plet) of SU(6) plays a central role [25]. This 
35-plet (35 = 1 + 24 + 5 + 5, in terms of 517(5) subgroup) has a VEV ~ Diag(l, 1,1,1, -2, -2). 
Therefore, in the context of F-like scheme it can provide a Clebsch factor k = —2, or, in T'-like 
scheme, k = 4 Q . 

In fact, one can consider more general cases of SU(N) group. Such a theory could contain 
adjoint Higgses with different structures which break SU(N) in various possible channels. By 
imposing the traceless condition on the SU(N) adjoint VEV, it is straightforward to obtain 
the relation between the two eigenvalues x and y in (PH ) and thus the relative Clebsch factors 
projected out in the context of the F-like and T'-like schemes. 

I. *i: SU(N) -» SU{5) x SU(N - 5) x [7(1), this Higgs contains only the SU{5) singlet VEV 
and thus can produce only k = 1. 

II. $ 2 : SU(N) -» SU(N - 2) x SU{2) W x [7(1), displaying the weak isospin SU{2) W . Hence 
($2) = Diag • (1, 1, — ^y^, — ^2~^) which implemented in the F-like scheme, generates the 
Clebsches k = -^f 2 , i.e. fc = -2, -5/2, -3, -7/2, • • • respectively for N = 6,7,8,9,- ■-. In 
the context of T'-like schemes the same VEV provides k = 2 ffi~ 4 2 ) , i.e. fe = 4, 10/3, 3, 14/5, 

III. $3: 5i7(iV) -» SU(3) c xSU(N-3)xU(l), exhibiting the colour SU(3) C . Therefore, ($3) = 
Diag-(1, 1,1, — 7v^3 ' ' ' ' ' — Tv^)' H ence 5 this Higgs in the context of F scheme gives k = —7^5, 
i.e. k = —1, —3/4, —3/5, —1/2, • • •, while in T'-scheme k = q^w, i.e. k = 00, —6, —3, —2, • • • 
(k = 00 should be understood as 1/k = 0, i.e. ^4 rf = while A e is non-zero). 

IV. There can be also adjoints with special VEV directions, having vanishing eigenvalues 
towards the SU(2) W components, say (<3?o) = Di a g"(T 1, 1,0,0, —-^5, • • • , — j^s)- It is natural 



to use such adjoints for the doublet-triplet splitting [26]. Obviously, this would give k = 
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in either F- or T'-like schemes. Analogously, there can be adjoints with vanishing SU(3) C 
components, which would lead to k = oo in F-like scheme and k = 2 in T'-like scheme. 

Therefore, in the next, while searching the realistic ansatze, we shall scan all possible Clebsch 
values for k n or l/k n among 0, ±1,±2, ±3 plus the above A-dependent expressions (for all 
possible A).0 

It should be remarked that within this wide variety, we find only three acceptable solutions 
which characterise three different ansatze. We anticipate them here: 

ansatz A k\ = k<i = 2 , k% = 1 , (59) 

ansatz B k\ = k% = —3 , ki = (60) 

ansatz C k\ = —3 , ki = , k^ = 2 (61) 

We shall come back to this point in Sect. 4.3. 



3 The 50(10) x SU(S) H model 

In the context of SO(W) model all fermion representations (|22|) are unified into 16-plets tpi, 
where i = 1,2,3 is a family (SU(3)h) index. In terms of the SU(5) subgroup, the fermion 
content is represented as tpi = (5 + 10 + l)j. So, in addition to the superfields 5j and 10j (p2|). 
tpi contain also the SU(5) singlets lj, usually referred to as RH neutrinos. 

The 5O(10) symmetry can be spontaneously broken to SU(5) by the Higgs supermultiplets 
C ~ 16, C ~ 16, having the VEVs towards the singlet components. The SU(5) symmetry can 
be further reduced down to SU(3) x SU{2) x U(l) by Higgses in tensor representations 54 and 
45. As for the MSSM Higgs doublets H u d needed for breaking the electroweak symmetry and 
generating the fermion masses, they are contained in a single superfield <p ~ 10.0 It is well 
known that the S'O(IO) framework also offers the possibility to work out the basic problem of 
the doublet-triplet splitting which has to be unavoidably adressed in any realistic grand unified 
theory. In particular, one can implement the so called missing VEV mechanism [^7]] which can 
be justified by additional symmetry reasons (see |28|, |29|] and references therein). It allows to 
achieve the mass splitting between the components of the Higgs 10-plet so that the doublets 
H U ( i remain massless while their colour triplet partners get masses of the order of the GUT 
scale, without unnatural fine tuning of the Higgs superpotential parameters. 

Since the 10-plet <fi can couple only to the symmetric combination of the fermion 16- 
plets, the higher-order operators inducing the fermion masses can only involve the (anti)sextet 
horizontal Higgs S (cfr. the SU(5) case ( p4[ ) where also triplet Higgses A can contribute). 
Hence, in the context of the S'O(IO) x SU(3)h model the operators (|J|) are represented as: 

O ~ ^Wi</>, (62) 

where, as in Sect. 2, the flavour scale M can be of the order of the GUT scale Mq ~ 10 16 
GeV, and Yt ~ 1 implies that £3 ~ M. One can still define S = M^ 1 (S). Thus, in the 

15 In the context of SU(N) models containing SU(5) x SU(3)h, (N > 8), the large value of the horizontal 
Clebsches b, b' could be similarly obtained if also the / + S ~ (1 + 8)-like combinations emerge from the VEV 
of the "big" adjoint breaking SU(N) down to SU(N - 1) x U(l), i.e. b = N - 1. 

16 We remind the content of SO(W) multiplets with respect to the SU(5) subgroup: 16 = 5 + 10 + 1, 
16 = 5 + 10 + 1, 54 = 24 + 15 + 15, 45 = 1 + 24 + 10 + 10 and 10 = 5 + 5. 
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low-energy limit this operator provides equal Yukawa matrices for all fermions, including that 
for the neutrino Dirac mass terms: Y u , Y^ e , Y^ 3 = S. In addition, the operator 

qij 

Or ~ w ^jC 2 . (63) 

involving the Higgs 16-plet with a VEV Vc ~ M induces the Majorana mass terms for the RH 
neutrinos: cx S. Thus, the effective operator for the neutrino masses emerges from the 
familiar see-saw mechanism, and it has a form (|5|) with Y„ = r/Y u , where r\ ~ (M/Vq) 2 - 

Therefore, in such a straightforward case we obtain a highly unrealistic pattern for the 
fermion masses and vanishing mixing angles in both the quark and lepton sectors. Recalling 
the discussion in Sect. 2.1, the way out from this situation is to invoke 50(10) and SU(5) 
breaking Higgses (45- and 54-plets) in the effective operators for fermion masses, as well as the 
SU(3)h octet Higgs for breaking the horizontal symmetry. 

In this section we present a consistent 50(10) x SU(3)h model to motivate the desired 
Yukawa textures (|20|). By explicitly specifying the extra symmetries, generically denoted by 
Q in the previous section, we will show, for example, how the Clebsch structure of the ansatz 
B may be obtained. Interestingly, the same symmetry reasons can be utilized for solving the 
problem of the doublet-triplet splitting via the missing VEV mechanism. 

Let us describe the ingredients of the model. For the sake of simplicity, we assume that 
all 'vertical' Higgses are singlets of SU(3)h and vice versa, all 'horizontal' Higgses are singlets 
of 50(10). In particular, the 'vertical' sector of 5O(10) Higgses includes chiral superfields 
in the representations 54, 45, 16 and 16, needed for the 5O(10) symmetry breaking down to 
SU(3) x SU{2) x U(l), and 10-plets for the electroweak symmetry breaking and the fermion 
mass generation. The "horizontal" sector of SU(3)h Higgses includes the triplet, sextet and 
octet representations. Namely, as in the previous section, we introduce one anti-sextet 5 ~ 6, 
three triplets ^-1,2,3 ~ 3 (plus the conjugated "spectator" superfields 5 ~ 6 and ^1,2,3 ~ 3), 
and the SU(3)h adjoint £ ~ 8. In addition, there can be some fields which are singlets under 
both 50(10) and SU(3) H - 

Besides the matter superfields tpi ~ (16, 3) we introduce a number of vector-like representa- 
tions containing the fragments like (^8|), ( p9|) and (|37|) which are needed for the fermion mass 
generation. Therefore, these should be matter superfields in 50(10) representations 16 + 16, 
10 or singlets, to be specified below. These states mediate the see-saw like diagrams and hence 
the quark and lepton Yukawa structures emerge after integrating them out. 

We invoke the additional symmetry in the form Q = U(1)a x Zq x 1Z, where U{1)a is 
an anomalous gauge symmetry []30| implemented in the same spirit as in ref. p9| , 1Z is a 
discrete symmetry (as in Sect. 2.2) under which all superfields in the theory as well as the 
superpotential change the sign, and Zq is a discrete symmetry acting on the superfields as 
<p n -» ^„exp(|i0„). 

The superfield content of the model with respect to 5O(10) x SU(3)h is given below, and the 
superfield charge assignment [Q, 0] with respect to U(1)a and Zq symmetries, respectively, is 
shown in subscript. Namely, we assume that the Higgs sector contains the following superfields: 

(i) the 5O(10) representations 

e~(54,l) [0)0 ] ) $~(45,l) [0i3] , n~(45,l) [0 , 2] , O' ~ (45, 1) [0)1] , 
C~(16,l) [ _ 1/2 , 1] , O- (16,1) [1A3] , ^~(10,1) [ _ 1 , ], (10,l)[i )3 ], 1 1 
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(ii) the SU(3)h representations 



S ~ (! 5 6 )[o,o]> Ax~ (l,S)[ 0t i], Ap. ~ (l,3)[o_ 2 ] ; A 3 ~ (l,3)[ 0il] , 

5 ~ (!> 6 )[0,0] ; ^1 ~ (1> 3 )[0,3]> ^2 ~ (l,3)[ 0) o], A 3 ~ (1,3)[ 0)3 ], 



(65) 



and E ~ (1,8)[_ 2( q], 

(iii) and the gauge singlets 

~ (M)[o,o]> ^ ~ (M)[0,2], Z~(l,l) [0 _ 2] ^ ~ (1, 1 )[-2,0]- (66) 

It is convenient to split the complete Higgs superpotential Wkiggs consistent with the above 
symmetries, into the following parts: 

Wi = a 2 y + y 3 + ye 2 + e 3 + (y + e)$ 2 + yzz + yss 
+z 3 + z 3 + zn 2 + zee + occ + zn' 2 + 

+S 3 + S 3 + ZA n A n + AxA 2 A 3 + A x AiA^ 
W 2 = W + I^ + ^C 2 , (67) 

(all SO(10) and SU(3)h index contractions are left understood). Here order 1 coupling con- 
stants are implied in the trilinear terms, and A is a mass parameter of the order of the GUT 
scale, Mq ~ 10 16 GeV.Q 

An important role is played by the D-term of the anomalous U(1)a symmetry 

TrO 

D A =M 2 + Y J Qn\Vn\\ M2 = T^ M * tT (68) 

where the sume runs over all scalar components <pk in the theory and Q are their U{\)a charges. 
Therefore, the spontaneous breaking scale of the U(1)a symmetry is naturally related to the 
string scale M str ~ 10 17-18 GeV with a coefficient determined by the charge content of the 
superfields in the game. 



Notice that A is the unique mass parameter contained in the superpotential (67), while 
another mass scale in the theory comes from the Fayet-Iliopulos term Ai 2 . Therefore, the 
VEV magnitudes of all fields into the game should be determined by these two mass scales, 
which in addition can be both of the order of the GUT scale. In particular, the Higgs fields 
participating in the superpotential W% can get ~ A VEVs violating the 50(10) and SU(3)h 
symmetries. In fact, the linear term A 2 Y induces the VEVs of the singlets Y, Z, Z which in turn 
play the role of mass terms for the superfields in non-trivial representations like f2, f2', <£, O, 
etc. Thus, with an accuracy of 0(1) coupling constants in the superpotential terms (|67[), the 
magnitude of all non-zero VEVs is ~ Mq. As for the fields I and E, they can get VEV of 
the order M. from the anomalous D-term (|6§|). We assume that the trace TrQ of the U(1)a 
charges over all superfields is positive: namely, the superfields presented in (^34|) and (^6|) plus 
the octet E contribute as TrQ ~ 10 2 . (The theory may include also other superfields with 
charges arranged so that the Green-Schwarz cancellation mechanism is at work, but this does 
not affect the orders of magnitude in our consideration). Therefore, Ai can naturally be about 
the GUT scale Mq or a bit larger. One has also to remark that the anomalous term fixes 

17 As far as R-parity is concerned, generally it is not automatic in 5*0(10) models involving the Higgs 16-plets. 
However, in the context of our model one can easily verify that all dangerous R-violating terms are suppressed 
by the Q symmetry. 
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only the overall VEV combination (I) 2 + (S) 2 , but not each term separately. The VEV values 
of I and S will be fixed by the soft supersymmetry breaking terms, as well as orientation of 
(E) with respect to the VEVs of other horizontal scalars, S and j4i,2,3) in the spirit discussed 
in sect. 2.2. This means that the horizontal fields should contain some light fragments, with 
masses order TeV, which, on the other hand, do not lead to phenomenological problems (see 
e.g. discussion in ref. [^]). The only important point is that, due to the simple structure of the 
SU(3)h group, the VEV of S in either case would acquire the simplest possible configuration 
(S)ocDiag(l,l,-2). 

Remarkably, the terms contained in the superpotential W\ provide the pattern needed to 
break SO(W) to the standard model group. In particular, there exists a solution for the 16-plet 
C, C VEVs which breaks £0(10) down to SU(5): (C) = (C) = V c ■ |+, +,+,+,+) (in terms 
of the corresponding Cartan sub-algebra generators). The condition (C) = {C) follows from 
the vanishing of the SO(10) D-terms. The 54 and 45 Higgses are needed to subsequently break 
£J7(5) down to SU(3) x SU(2) x U(l). The necessary VEVs are given as: 



(0) 


= V e 


• Diag(l,l, 1,-3/2, -3/2) ® <x 


<*) 


= Vq> 


Diag(l,l,l,0,0) ® a, 


<n> 


= V n 


Diag(l,l,l,l,l) ®a, 


<n'} 


= 0, 





(69) 
where 



One can easily verify that these VEVs are consistent with a super symmetric vacuum and their 
configurations are not affected by some unknown factors. In particular, in the 54-plet G only 
its 24 fragment (in terms of SU(5) subgroup) acquires VEV and, noticeably, into the only 
phenomenologically allowed direction which induces the 50(10) symmetry breaking down to 
the Pati-Salam subgroup: 50(10) -> ££7(4) x £17(2) x SU{2)'. The coupling of the 45-plet 
<3? to in the superpotential forces the VEV of the former to lie in the B — L direction, as 
a result of a combination of its SU(5) 1 and 24 fragments, with exact zeros on the last two 
components, as it is demanded by the missing VEV mechanism |27j , pgfl.P] As for the other 
45-plet f2, it is coupled to the 16-plets C, C, but does not couple to the 54-plet 0. Because of 
this, its VEV is directed entirely towards the £f/(5)-singlet direction, with a vanishing VEV 
towards the 24 fragment .0 And finally, the third 45-plet f2' does not couple neither to nor 
to C, O, and thus it should have a vanishing VEV. 

It is also important that all non-singlet (under the standard model group) fragments in the 
£O(10)-breaking superfields get masses of the order of Mq. In such a case, no light fragment 
(with mass <C Mq) is left behind which could affect the gauge coupling unification at the scale 



Mq. The presence of the term <&Q£l! is crucial for achieving this result [28|, since otherwise 



in its absence, the theory would have flat directions towards non-singlet fragments which are 

18 As we will see below in the analysis of the Yukawa sector, such a VEV structure of $ is also important for 
obtaining the Clebsch factors needed for the ansatz B, namely fci,3 = —3 and fe2 = 0. 

19 This is also relevant for the fermion Clebsch structures, because a S l l7(5)-singlet VEV - namely, the combi- 
nation (Z + Q) - participating in the fermion mass generation (see below, diagrams in Fig. ^) does not induce 
uncontrollable corrections to the Clebsch factors. 
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related to extra global symmetry in the superpotential. Remarkably, this term does not affect 
the VEV pattern of 45-plets while lifting the dangerous flat directions. 

Finally, the term W\ in ( |6T|) is also responsible for the horizontal SU(3)h breaking via the 
VEVs of S and -Ai,2,3 fields (see ( |26| ) and (p7|)) and its features are elaborated in Appendix 
A. Notice, that even if all coupling constants are taken real in the superpotential, non-trivial 
CP-violating phase in these VEVs could appear in case of solutions with different phases for 
the singlet VEVs (Y) and (Z), in the spirit discussed in Appendix B. 

Let us address now the doublet-triplet splitting problem. As well-known, the VEV of the 
54-plet guarantees the VEV of the 45-plet $ to be towards the B — L direction, ($») oc 
Diag(l, 1, 1, 0, 0) x a. On the other hand, the superpotential Wi contains the necessary terms 
to incorporate the missing VEV mechanism for the doublet-triplet splitting. In addition, they 
also induce the mixing between the 5 fragments in 4> and C through the large SU(5) conserving 
VEV Vc = (C). All this can be more easily represented in terms of the following mass matrix: 



(71) 



where M ~ (/) and Mc ~ (Z + Q) is the mass of the 5-plet fragments contained in C, C. 
Owing to the form of the VEV ($) in (|6^), non-zero entries are induced only between the 
triplet fragments in (j> and 4>' while the corresponding entries for the doublet fragments vanish. 
Therefore, all the triplets get masses of the order of Mg, while two combinations of the doublet 
components, to be identified with the MSSM Higgs doublets H u ^, remain light. The doublet H u 
comes entirely from H = 5(4>), and so no portion of it is contained in 5(C) or 5(c/>'). As regards 
Hd, it emerges from the combination of 5 fragments in and C. Namely, H = Cuj5(C) — (0) , 
with tana; ~ Mc/Vc (see also PH). 

Let us discuss now the fermion mass generation mechanism. Apart from the already familiar 
tp ~ (16, 3)[_i/2,o] ; let us introduce the following vector-like fermions: 

(I) 16-plets ' 

*~ (16, 3) [3/2,0], ~(16,3) [1/2i0] ; *' ~ (16, 3) [3/2)3] , ~ (16, 3) [1/2)3] , (72) 

(II) 10-plets 

t~(10,3)[_ lj3]j *~ (10, %,_!], ti~(10,3) [lf2] , ti~(10,3) hl . 2] , 

*2~(10 J 3) [ _ 1 ,o ] , i 2 ~(10,3) [liQ] 

(III) and singlets 

iV~(l,3) [0j3] , F~(l,3) [0i3 ]. (74) 

Let us start by discussing the role played by the additional fermion 10-plets ti = (5 + 5)j, 
? = (5 + 5) 1 . All the symmetries of the model allow in particular the following terms: 

Cipi+Ztl (75) 

Correspondingly, the mass matrix connecting the matter states reads as: 



(73) 





5(*i) 


5(f) 


5(&) 




V c 


5(ti) 


(? 


(Z) 


m 


V (z) 






(76) 
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Figure 3: The see-saw diagrams giving rise to the effective higher order operators for the symmetric 
contributions to the charged fermion Yukawa matrices and to the neutrino Dirac mass term Y^ 3 - 
diagram (a), and to the Majorana mass matrix Mr of the RH neutrinos 1 contained in ip — diagram 
(b). In SU(5) language the relevant fragments are ip — 5 + 10+1, = 5 + 5, and the singlet component 
of C with VEVVb. 



which shows that the light states 5j are composed by a combination of 5 fragments in ipi and 
U i.e. hi = — S0b(ip)i + c$5(t)i, where the angle 9 is defined as tan# = (Z)/Vc, while the 
orthogonal combinations 5^ = cgb(tp)i + Sg5(t)i and 5(t) 1 are heavy states. 

Now we are in position to fix the fermion mass pattern which emerges in our theory. The 
diagonal entries in the Yukawa matrices are generated by the following superpotential terms: 



Wm = ip^S + *(/ + £)* + ^Sip 
W N = ipCN + NYN + SN 2 , (77) 

which through the see-saw like diagrams shown in Fig. 0, give rise to the effective operators 



like (32) and (p3[). For an easy identification of the contributions to the Yukawa matrices, the 



terms in and Wjy can be rewritten according to their SU (5) decomposition: 

Wis, = 1O(^)5(0)1O(«) + 5(^)5($10(tf ) + ID(*)510(^) + 10(*)(J + E)I0(tf) + 

10(^)5(^)5(*) + 5(V)S5(ip) + 5(tf)(J + S)5(¥) + 

5(^)5(0)1(*) + + )(/ + S)l(¥), 

Wiv = l(tp)l(C)N + NYN + SN 2 . (78) 

Therefore, as already discussed in Sect. 2.1, the VEV combination (J+S) provides the SU(3)h 
breaking and leads to the SU(2)h invariant form of the heavy fermion mass matrix M$ = Mb', 
where b' = Diag(l, 1,6') and M ~ Ai. The relevant Yukawa terms for the light fermions emerge 
after integrating out the heavy states. In particular, for the Yukawa coupling matrix of up 
quarks (which are entirely contained in 10 (ip) fragment) to the Higgs doublet H u C 5(<p) we 
obtain Y n = M^S 1 ) = b /_1 S, where S = M^iS) and it is diagonal once the VEV (S) is 
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Figure 4: Diagrams giving rise to the effective higher-order operators for the non-diagonal contributions 
to Y e: d. The relevant contributions of these diagrams involve Hd C 5(C) and 5, C 5(tj) states. 



chosen as (26). On the other hand, recalling that 5(i(i) contains the light 5-plet with a weight 



sg and the Higgs 5((f>) contains with a weight s^, the terms (ff^) provide also the diagonal 
contributions to the down quark and charged lepton Yukawa matrices Y^ s = sgs UJ h f ~ 1 S. 
From the decomposition (|7^) it is apparent that the (heavy fermion) exchange contributions 
to Y u are of the T -scheme type (cfr. ( |35| ) ) and those to Y^ e s are of both T - and F -scheme 
type (cfr. also (|38|)). 



Finally, for the Dirac Yukawa matrix of neutrinos we obtain Y^* = seb /_1 S. On the other 
hand, by considering the terms in Wn, we see that the VEV (Y) induces the SU(3)h invariant 
masses between the states N and N, Mjv = Mn • 1, (1 = Diag(l,l,l)) and ~ Mq, 
and so the RH neutrinos l(ip) acquire the Majorana masses M R = M N 2 Vg(S) = kMS, where 
k ~ (Vc/Mjv) 2 ~ 1. Thus, as long as both Dirac and Majorana terms are present for neutrinos, 
the operator (|J) emerges via the standard see-saw mechanism: 



Y s 2 
M L ~ Y " M * Y " ~ kM 



YjfM^Yj? = -^jb'- 2 S (79) 



Therefore, both Y n and Y v are diagonal and we have: 

Y„ = rjb'- l Y u (80) 

where r\ = Sq/k. 

Let us discuss now the origin of the off-diagonal entries in Yd e , recalling that they have 
already got the diagonal contributions tx Y u . We also remember that the off-diagonal entries 
should contain non-trivial Clebsches between quarks and leptons. 

To this purpose, we recall that the light fermion 5-plet is also contained in 5(t) with a 
weight cq, while the doublet is contained in the fragment 5(C) with a weight c w . Then, 
the symmetry of the theory prescribes other superpotential couplings for the superfields t and 
C, involving other (than ^, \&) heavy fermions, so that the exchange of the latter induce the 
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off-diagonal contributions to Y<2 >e via the VEVs of the triplet horizontal Higgses -Ai,2,3- The 
relevant diagrams^ are shown in Fig. ||. The diagram (a) provides the entries (2,3) and (1,2) 
as only the SU (3)# triplets A\ and A3 are allowed to have a direct coupling to the matter 
10-plets ti. In this case, the coupling of the 45-plet $ (with its VEV into the B — L direction) to 
the 16-plets (ifi and ^ ) gives rise to the Clebsch coefficients ki t 3 = —3. We should also observe 
that these off-diagonal contributions are effectively induced by the exchange of the SU(5) 10+ 
10 fragments contained in the vector-like states + , according to the T' -scheme (cfr. (flO|))- 

The triplet A2 is instead involved in the diagram (b) which therefore contributes to the 
(1,3) entries. Because of the requirement of Q invariance, this triplet cannot couple to the 
states ti but it couples to the heavy states ti which are mixed to the former through the VEV 
of the 45-plet This occurs through the couplings on the left of the A<i incoming- line in the 
diagram (b), while the part on the right side contains exactly the same couplings as that on the 
right of the A\ % line in diagram (a). However, the coupling of the 45-plet $ with the fermion 
10-plets, t&fa, gives vanishing contribution to the lepton doublets contained in 5(i) (similarly 
to what happens in the missing VEV mechanism for the doublet-triplet splitting) and therefore 
we obtain k<i = 0. In either case of diagram (a) or (b), the SU(3)h breaking Clebsch b emerges 
from the same combination I + S (i.e. from the matrix M^, ~ Mb^/ = M ■ diag(l, 1, 6*')) so 
that the 13 and 23 entries have the same asymmetry. 

Thus, the final form of the down quark and charged lepton Yukawa matrices reads as: 

Y d = pY° + b" 1 A d , Y e T = pY° + b^Ae, (81) 

where the anti-symmetric matrices A^ e have the form (j50|) with = —3 and &2 = 0, 
b = diag(l, 1, b) where b = fe^, 1 . In this way we have reproduced the Clebsch factors demanded 
by the ansatz B (|60|). The coefficient p is identified with s^sg. It may be smaller than 1 and 
thus moderate or even small values of tan/3 can naturally be accommodated. 

In conclusion, we have presented a complete 50(10) x SU(3)h model in which the symmetry 
breaking VEV pattern as well as the form of the superpotential couplings needed to obtain 
the desired predictive mass texture for the fermions (specifically, the ansatz B presented in 
Sect. 2) are fully motivated by the additional Q = U{\)a x Z§ x 1Z symmetry. Interestingly, 
the missing VEV mechanism for the doublet-triplet splitting can be motivated on the basis of 
the same symmetry Q. The other Clebsch structures demanded by the ansatze A or C can be 
reproduced along similar lines by rearranging the charge assignments of the superfields in the 
game with respect to the Zq symmetry, however we will not show the corresponding models 
here. 

4 Phenomenological analysis 

This part is devoted to the phenomenological study of the ansatze suggested. First, we provide 
the corresponding expressions of the physical quantities - quark masses as well as mixing angles. 
Second, the criteria of the fit are presented. Finally the results are shown and discussed. 



For the sake of brevity we do not write down the relevant superpotential terms responsible for the generation 
of the off-diagonal Yukawa couplings: they can be read out directly from the diagrams in Fig. W. 
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4.1 Physical observables from the theoretical model 

We now deal with the explicit form for the Yukawa matrices ~Y e ,d in (|5ll). In general the 
coupling constants and the VEVs present in the theory are complex. However, all the phases 
except three can be pulled out by field redefinitions and the Yukawa matrices can be presented 
in this forrmPI 

/ >« 

Y c I , Y v = n Y c I , (82) 




V 




(83) 



/ Y ut De l( Ae i(J 
Yd = -Ae™ YaDe* 
V \C \B 

Here D = pY t , Y ut = Y u /Y t and Y ct = Y c /Y t . Notice that in Y^, b f — b has been taken 
for concreteness. As compared to the expressions in ( |5l|) the off-diagonal terms have been 
redefined as Af = B, A^ = C and Af = A and understood to be real.0 

The following observations can be made. The 11 entry in the matrices Y e ^ is negligible 
and in the following we set it to zero. Indeed, the ratio Y u /Yt is two orders of magnitude less 
than Y e /Y T and Y^jY^ (cfr. Table [j]) and thus its contribution in Y^ e is irrelevant. 

As regards the 22 entry, this also leads to small corrections - the ratio Y c /Yf is one order of 
magnitude less than Y s /Yb (cfr. Table [l]). However, we find that these corrections are relevant 
on comparing with the present precision in the CKM angles. So this contribution in Yd e cannot 
be neglected. 

In general case the phases in ( p3| ) are arbitrary. However, they can have specific values 
in the context of the spontaneous CP violation. Namely, one can assume that the original 
theory has an exact CP invariance, i.e. all couplings in the Yukawa and Higgs sectors can be 
simultanoulsy made real by phase redefinitions of the superfields. However, CP invariance can 
be spontaneously broken by non-vanishing phases of some VEVs in the theory. In particular, 
as we show in Appendix B, the triplets A n can have complex phases. 

For example, one can choose a basis where the VEVs of the sextet S, 5i,2,3 are real. Thus, 
in this basis all three diagonal entries of Y^ in ( |5l| ) are real. For definiteness, let us take them 
also to be positive. On the other hand, the off-diagonal entries keep on the phases originated 
from the VEVs of A n . As we show in Appendix B, these VEVs can exhibit a spontaneous 
CP violation when the product ^1^42^3 is imaginary. In particular, if all entries ^1,2,3 are 
imaginary, the matrix Y^ in ([5l]) , once recasted in the form (|83|) , implies the phase assignment 
£ = vr, a = tt/2. In another possible case, when A12 are real but ^3 is imaginary, one has 
instead £ = and a = tt/2. In principle, the spontaneous CP breaking could provide some 
other interesting phase assignments such as £ and a have some "inverse 7r-fold" values like 
2ir/3, 7r/4, etc. In the following, to be more general, we keep the phases £ and a arbitrary, 
however we shall pay special attention to those specific phase values which might be obtained 
in the context of the spontaneous CP violation. 

21 As far as Y u and Y„ are diagonal, these phase transformations are not relevant for the form of the mixing 
matrices V q ,i and we do not show them. 

22 In general the relative "lepton versus quark" Clebsches k n — A^/A 1 ^ can be complex, however in the 
following they are taken real. 
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Let us now proceed our analysis. The (1,3) and (3,1) entries in Y d ^ e can be simultaneously 
rotated away by orthogonal and phase transformations in the 12-flavour space: 



P 



^Y d O d 12 P 




P 



0\ 2 P. = Y e , 

/ ? d{e) ~d(e) 



n d ( e ) 

u 12 



12 



V 



12 



d(e) d(e) 
s 12 c 12 





(84) 



where c 



d(e) 
12 



COS 



7 12 



etc. 



and it is taken tan^^ 2 



C/B and tan#f 2 



^tan0 d 



fei 



12- 



Notice 



that these rotations do not affect the (1,2) and (2,1) elements as the Yukawa matrices are 
antisymmetric in the 12 block. On the contrary the 23 block gets modified and by further 
re-defining everywhere the parameter B as B = ■—[-, the matrices become: 



/ -.4 
-A Y ct De^ B 
{ \B D 



( -k 3 A 

-k 3 A Y a De< fhB 
\ k\B D 



(85) 



where k\ =k\J (cf 2 ) 2 + (sf^^/^i) 2 Finally, these matrices Y e and Y^ can be diagonalized by 
bi-unitary transformations as in eq. (|6|) with the matrices U ejd and U' e d understood as: 



u = u 32 u 13 u l2 = 

( 1 

C 2 3 
V -s 2 3e~^ 23 




-si 3 e 


C13 



-i<t>l3 



\ 



C12 





-i<t>i2 



sue 
C12 




(86) 





1 



Therefore, the CKM matrix is V q = Of 2 P + U d and is parameterised as 



V„ 



( cf 2 e" i 

V 



v ud ~\~ ^i2 v cd ^12 e ^ v us T S \2 L < - 
■xd „—ia„, ~d „, ~d „—ia„ 



'12 e v ud C\2 v cs S 12 f 
V t d V ts 



p,d „—tu 

c 12 e 

c l2 v cb 



Vub + Si 2 V c b 

- s 12 e la v u ij 

Vtb 



(87) 



where by v ud , v us etc. with lower v we mean the matrix elements of the quark mixing U d 
as defined in eq. (R) and parameterised as in (|8|). Analogously, the effective leptonic mixing 
matrix is V t = U^PLOfJ: 



( 



-c\ 2 e la v el + s\ 2 v e2 c e 12 v e 2 + s\ 2 e la v e i v e3 
-c\ 2 e lcJ v^i + s\ 2 v^2 c\ 2 v^2 + s^e^Vfj,! v^ 3 
\ -c\ 2 e m v Tl + ~s\ 2 v r2 c\ 2 v r2 + s\ 2 e la v T i v r3 



where with lower v e i,v e2 etc. are denoted the matrix elements of the leptonic mixing matrix 
U\ . These forms of quark and leptonic mixing matrices are convenient since the contributions 
v U s,v ud and v e \, v e2 etc. from the Yukawa matrices Y e ^ d of (|85|), are separated from that arising 
from the initial 12-rotation 0^ 2 6 . Had the matrices Y £jd the form in ([i~3| ) (i.e. vanishing 22 
entry) the above transformations U ed would be pure rotations. Then the only source of CP 
violation would reside in the residual phase a in P± driven by the (1,3), (3,1) entries. It 
is worth noticing that the CP violating phase appearing in the quark and leptonic mixing 
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matrices differs by it (cfr. (§))• This feature comes from the antisymmetry of the 12 block, 
while it does not depend on the initial antisymmetry in the 23 and 13 blocks. 

The GUT scale Yukawa eigenvalues Y u>Ci t, Yd >S) b, Y e>ll>T , ^1,2,3, are linked to the physical 
fermion masses through the renormalization group equations (RGE). For moderate values of 
tan/3 = v u /vd, one obtains at one-loop^] (see e.g. p2|): 



m u — Y U R U R U ^ T] u Bfv u , 
m c = Y c R u R^ M Vc Bfv u , 
mt = Y t R u R s u M Bfv u , 



m d = Y d R d R s d M Vd v d , 

m s = Y s R d R% M r] s v d , m v 

m b = Y b R d R% M r] b B t v d , m T 



Y e R e R s e M v d 
: Y^R e R s e M v d 
Y T R e R^ M v d , 



(89) 



and for the neutrino masses 



mi, 2 ,3 



M v v 



SM R 6 2 
B~ t v 2 



(90) 



where the factors R u ,d,e,i/ and R^Jde v account for the gauge-coupling induced running from 
the GUT scale M G ~ 10 16 GeV to the SUSY breaking scale M s ^ M t and from M s to the 
electroweak scale Mz, respectively. The factors r]f encapsulate the QCD+QED running from 
Ms down to mj for f = b,c (or to fx = 1 GeV for the light quarks / = u,d,s). Namely, for 
a s (M z ) = 0.119 ± 0.004 we have 



E> DO 



SM 



Vb 
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K d K d 
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+0.13 
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H e H e 
Vu,d,s 



1.50, 



2.38 



-0.24 
-0.19- 



R V R V 



SM 



1.15, 



(91) 



The factor Bt includes the running induced by the large top quark Yukawa constant (Yt ~ 
1).0 For Yt varying from the lower limit Yt = 0.5, imposed by the top pole-mass, to the 
perturbativity limit Yt ~ 3, the function Bt decreases from 0.9 to 0.7. Regarding the CKM 
elements, their physical values are related to the corresponding GUT-scale quantity (labelled 
by the superscript G) as follows 0: 



V, 



us(d) 



V us(d) 

vl.^B: 



c(u)b — v c{u)b D t 



cs(d) 



td(s) 



V G 
V cs(d) ' 

Vtd{s) B t 1 



Vt 



Lb 



tb 1 



(92) 



Using the procedure outlined in detail in [ 13 1 , the following relations (valid at the GUT scale) 
are derived from the matrices (pq): 



D = Y T 



1 



1/2 



Y b 



1 



(b + b-^gg 
F~d 



l + cf 2 )2 Y c t{b + b~ l ) C osi 



1/2 



(93) 
(94) 



and 



A 2 D 



YY Y 
1 e 1 11 1 t 



HY d Y s Y b 



(95) 



The notation is the standard one: for the heavy quarks t, b, c, mt,b,c are their running masses respectively 
at /i = mt,b,c, while for the light quarks m u ,d. are given at fi = 1 GeV. 

24 In the RG running of neutrino masses (|9(]|), the factor Bi is function of Yt/ sin f3. 

25 We shall not take into account the analogous RGEs for the neutrino mixings, since the experimental 
data still contain big error bars and such an improvement is not justified. Moreover, renormalization effects are 
mostly expected in case of strong mass degeneracy |34[ which is not our case. 
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j- 2 {Y^ - Y e )Y T f = b(Y s - Y d )Y h -± , 



Id 
'F d 



where we have defined 



1 + cf| 



Y T 



Y — Y 

1 n 1 e 



Y a cos £ 



I d =l+ c 12 



Y - Y, 



Y ct cos i 



(96) 



(97) 



In the following we directly substitute the Yukawa constant ratios with the corresponding 
mass ratios whenever the latter are RGE invariant, e.g. Y^/Y T = m^/m T , Y d /Y s = md/m s , 
etc. Now, by dividing the squared of the l.h.s and r.h.s in eq. (|93|) by the corresponding sides 
of eq. (p6h and inverting, we obtain the following implicit relation: 



where 



Z 



Y s -Y d 



F P 



IeF^ 2 



Trip 



I d F e Z 2 



0.059 



1 



1 



TO, 



I e F d k{ 2 
I d F e Zi ' 



, 1 x m u — m e 



TO T 



1/2 



Substituting the expression (^g) back into eq. ( p3[ ) we get: 

Z, 



Yb 
Y T 



(98) 



(99) 



(100) 



which shows the modification induced on the expression of the b-r Yukawa unification at the 
GUT scale by the asymmetry in the 23 block. Analogously, by utilizing the relation (100) in 
eq. (|96|), the latter can be rewritten as 



fcl IeF d 

Z I d F e 



(101) 



Y^-Y e 

while by dividing the squared of the r.h.s. and l.h.s. of ( |96[ ) by the corresponding sides of (^5| 
we have: 



m s | m d 
m d m s 



k 2 



hF d 



kfZ \I d F e 



+ — 

m P . to. 



204.7 



kl (I e F d ^ 2 



kfZ \I d F e 



(102) 



Compare the present expressions ( 38 , 101 , 1Q2| ) and ( |99"|) with the corresponding ones in eqs. 
(14,1^) and (O) obtained in the Fritzsch-like case (C = 0, Y ct ,Y ut — > 0, i.e. F e ^l e ^ d — > 1). 
From eqs. (liOCl ) and (|l0l| ) we can extract the following physical masses: 
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We now turn to the quark and lepton mixing. For the angles 6^ and the corresponding 
phases appearing in the unitary transformation (86) we find: 
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(106) 
(107) 



while the right rotation angles an d 023 are obtained from these expressions substituting 
respectively b — ► The phases appear strongly suppressed by the tiny y rf parameter: for 
Y ct ~ 3 • 10~ 3 and 6 ~ 10, < 0.01 . The 13 mixing angles are 
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One can easily verify that for the corresponding phases </> 13 
mations are expressed by the following angles: 
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Finally the U\2 transfor- 



tan20 



0.140 



12 



tan 26>f 2 



while the corresponding phases are 
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(113) 



We can easily estimate that (f>f 2 can be at most 11° and 25° with b > 1, yt = 3 x 10~ 3 and 
for ki = 2 and fci = —3 , respectively. Moreover, the larger is the asymmetry parameter b, the 
smaller the function Z becomes and therefore the smaller the phase <^>f 2 appears. In conclusion, 
the phase £ cannot provide by itself a sizeable source of CP violation since it is driven by the 
tiny ratio Y c t. That is the reason motivating us to introduce non-vanishing (1,3), (3,1) entries 
driving the extra phase a. Hence, as a matter of fact we may think of the matrices U e d as 
pure rotations, parameterised in the standard form (]7|) with the phase 5 = tt jOJ. 

In the analysis, we test the amount of CP-violation through the parameter ek describing 
the CP violation in K —* hit. More precisely, we use the more uncertain parameter Bk 
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Observables 


Values 


• m e [MeV] 


0.511 


• m M [MeV] 


105.7 


• m T [GeV] 


1.777 


• m u /m t 


(0.9 H- 3) x 10~ 5 


• m c /mt 


(8.0 ± 1.3) x 10- 3 


• Q 


22.7 ±0.8 


• \ £ k\ 


(2.280 ± 0.013) x 10" 3 


m s [MeV] 


155 ± 75 


m b [GeV] 


4.25 ± 0.15 


771- / 777 j 


21 ± 4 


TYl„, I TtlA 
" b U/ " "a 


0.2 ~r~ 0.7 


M t [GeV] 


173.8 ±5.2 


1^1 


0.0395 ± 0.0017 


\v us \ 


0.2196 ± 0.0023 


\v ub /v cb \ 


0.093 ± 0.014 


\Vtd\ 


0.0088 ± 0.0018 


\Vtd/V u \ 


< 0.22 


B K 


0.8 ±0.2 



Table 4: Physical quantities used in the phenomenological analysis. Those marked by • are used as 
input for the fits. The bounds on |Vtd| and |Vtd/Vf S | are inferred from the evaluation of 73° — B® and 
B® — Bg transitions in the standard model || and thus they are not literally valid in the context of 
supersymmetric model; even in the MSSM framework, with flavour-aligned soft terms, one can expect 
substantial supersymmetric contributions to these transitions. 



parameterising the deviation from the vacuum saturation limit in K° — K° transition which 



we 'extract' from the expression of £k as given in the standard model [35]: 
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n -1 



(114) 



and compare with its theoretical value which ranges from 0.6 to 1. Here fx = 161 MeV is 
the kaon decay constant, M K = 497.7 MeV is the K° mass, and AM K = 0.53 x 10 10 s _1 
is the experimental value of Kl — K$ mass difference. The mixing angles enter through the 
combinations t a = V* s V a d'- t u appears explicitly,^] whereas t c , t t are contained in the function 



T = rjitlSo(x c ) + mttSo(xt) + 2r] 3 t c t t S (x c , x t ) 



(115) 



where r/i ~ 1.38, 772 ~ 0.57, 773 ~ 0.47 are QCD correction factors and is one of the Inami- 

2 

771 

Lim functions (x c t = — I 2 -). In this way, we are testing whether the experimental value of 

' m w 

£k can be reproduced entirely by the standard model contribution associated with the box 



26 Note that in most of Vckm parameterizations, the factor t u , coming from the amplitude of the decay 
K° — > 2ty, is made real and thereby the factor £„/|tu I becomes 1. In our final parametrization of the CKM 
matrix this is not the case - see (p7|). 
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diagrams involving the charm, top quark and the W boson. It is known, that in the MSSM 
or more generally, in models with flavour-aligned soft terms like in []36| , the supersymmetric 
contributions are negligible. However, in general SUSY GUT context ek can receive substantial 



contributions from the super-partners at the weak scale [BT], 38]. In view of this, we shall not 



use the information from ek to constrain the parameter space. 



4.2 Model parameters versus physical observables 

We can now proceed by confronting the parameters with the physical observables described by 
the model. The number of the former amount to 14, namely Y u>Ci t, A, B, D, if 2 , b, a, £, ki, k 2 , 
k%,ML while the number of observables is 20: the up-quark masses Tfiuctt the down-quark 
masses m,d, s ,b, the charged-lepton masses m eiMjT , the neutrino masses 7711,2,3, four independent 
CKM parameters and other four parameters of the leptonic mixing matrix. As it will be 
discussed more extensively below, the three Clebsch coefficients ki,k 2 , k^ are treated as external 
parameters, fixed at some specific GUT-inspired value. 

The various physical observables deserve a different treatment in the analysis for at least 
two reasons. First, they are known with different degree of accuracy. Second, some of them 
have a minor impact on the determination of the 'flavour' parameters themselves. Therefore, 
the observables are classified and used as follows: 

1. m e , m^,m T : these masses are determined with high accuracy so they are used as inputs to 
extract the parameters A, B, D via eqs. ( |93"y96| ). Notice that in this way, also the explicit 
dependence on tan (3 of the fermion masses is automatically absorbed in the parameters 
A,B,D. 

2 mu me. th ese mass ratios enter only into the determination of and are used as 

inputsp 7 !- Notice, however, that Y u t,Y c t quite depend on Yt through the renormalization 
factor B t : e.g. Y ct = mc/mtB^r]^ 1 . As this ratio is tiny, the dependence on m c /m t is 
weak and hence our results are not sensitive to the corresponding errors on m c and nit- 
This justifies the fact that we use the central values of the ratio m c /mt shown in Table 
|j. On the other hand, the Yukawa constant Yt is left as free parameters. Indeed, it is 
not sensible to fix Yt by the top mass due to its infrared behaviour. It will turn out 
that (in almost every case) the preferred value for Yt is 0.5 which enhances the effect 
of Y C f. This value is indeed the smallest one compatible with the top mass for which 
m t = (160.9 ± 2.2) sin (3 GeV for a s (M z ) = 0.119 ± 0.004. Correspondingly, we have for 
the ratio Y ct = (2.7+^) ■ 10~ 3 . 

3. rrid, m s ,mb, \ V US \, \V c b\, \V u b/Ycb\'- these observables are fitted by a x 2 — like procedure (see 
Sect. 4.4) to constrain the remaining parameters they depend on, namely Yt , b, if 2 , c, £ 
once the Clebsch factors ki,k% are prescribed (see Sect. 4.3). 

4. M t ,m u /m^: the top pole mass M t = m t ■ (l + ^a s (M t ) + ^£a 2 s {M t )^ is determined 
as output once Y t is fixed by the fit in point 3., modulo sin/3. Though in our analysis 
we cannot determine univocally (and in a sense is not necessary) tan f3, the model we 
presented in Sect. 3 naturally yields tan/3 in the small or moderate regime. Using the 

27 Now for the sake of discussion we consider also Y u t but in practice, as already mentioned, we neglect it at 
all. 
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determination of m s /md from the fit, we extract the ratio m u /md from the combination 

n 

m s /m d 

Q= n 7= i y 2 - ( 116 ) 

VI - {m u /m d y 

5. Bk- the CP violation parameter is given as a prediction from the result of the fit per- 
formed in point 3. 

6- \V^\, \ V e 2\, \ V e 3\ & lepton CP-phase: these are also given as predictions. More precisely, 
for the sake of comparison with the experimental data, we trade these three elements 
for the corresponding 'oscillation' parameters: sin 2 29 l 23 = 4| V^3| 2 (l — (V^l 2 ), sin 2 26[ 2 = 
A\V e2 1 2 {l-\V e2 1 2 ) and sin 2 29 l 13 = 4|V e3 | 2 (l- 1 V e3 \ 2 ). Notice, however, that the predictions 
for both CP violation phase and the element \V e2 \ depend also on t\ 2 = ^t\ 2 i.e. on the 
Clebsch k 2 . 

7. mi, m 2 ,m 3 : the neutrino masses follow immediately as predictions m 2 /m 3 = bY c t, 
m\/m 2 = Y uc . However, the determination of any single mass does require the last 
parameter i.e. the mass scale Ml = —. 

In summary, we have that from the sets 1) and 2) five parameters are fixed - A, B, D, Y u t, Y C f. 
From the set 3) other five parameters are constrained by the fit - b, tf 2 , a, £, Yj once the two 
Clebsches k\, k% have been assigned. The determination of these twelve parameters allows us to 
make five clean predictions: nii/m 2 ,m 2 /m 3n Bk, IV^I, \ V e $\. Also M t is given as a prediction 
within the uncertainty of tan (3 expected to be in the range 1 — 10. From the result of the fit, 
we also gain the ratio m u jnid by using the information on the parameter Q (fixed at its central 
value) . 

Finally, the neutrino masses depend on Ml and both \ V e2 \ and the amount of CP violation 
in the lepton mixing matrix depend on the Clebsch coefficient k 2 . For definiteness, the flavour 
mass scale M and the lepton- violation mass scale Ml are set at the GUT scale Mq, i.e. M = 

Ml = Mq which as already discussed means that — — ~ 10 or n ~ 10 2 for Y t ~ 1, b ~ 10. 

rfY t 

In the next Section we shall discuss more extensively how we handle the Clebsch factors. 

Finally, the phenomenological analysis requires also the parameters associated with the 
standard model gauge sector: a s (Mz), a em , sin 2 9 W . For a em , sin 2 6 W we use the central values 
quoted in ||. As regards a s (Mz), the present world average is a s (Mz) = 0.119±0.002 Q. On 
the other hand the SU(5) unification of gauge couplings implies a larger value, a s (Mz) ~ 0.125. 
Therefore, as a compromise, we have taken a s (Mz) = 0.123. We remind that only nib and m s 
are sensitive to variations of a s (Mz) within the quoted error. The corresponding variations 
can be inferred from the gauge-renormalization factors given in (pl[). On the other hand, the 
a s (Mz) dependence of the running factor Bt is less than 1% and therefore all other observables 
are not sensitive to a s {Mz)- 



4.3 Clebsch prescription: three different ansatze 

First let us recall how the Clebsch coefficients enter into the determination of the physical 
quantities. The Yukawa eigenvalues Y e ^ jT and ld,s,6 depend only on k% and k% (see the analyt- 
ical expressions in Sect. 4.1)0 Similarly, the mixing angles v us ,v u d, ■ ■ ■ also depend only on 

28 The dependence on ki contained in k\ in eq. (EH) is mild for tf 2 < 0.1. 
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k\,ks- The coefficient k 2 affects only the lepton mixing matrix (|88[ ) via the initial 12 rotations 
t\ 2 = j^tf 2 . Thereby k 2 can only be fixed for example by the 12 lepton mixing angle V e2 
through the SN parameter sin 2 29[ 2 = 4|V e 2| 2 (l — IV^] 2 ). Notice, indeed, that both sin 2 29[ 3 
and sin 2 29 2 % do not depend on k 2 . 

We can guess very easily the needed range for such k± 23- Indeed k% t 3 are mainly fixed by 
the strange mass in ( |103D and the ratio m s /md in ( |102| ). From m s we infer that k\ is to lie 



within ! -j- 3 - for 6 > 8. Correspondingly, from m s /m^, ^3 is forced to be in the range 1 4- 1. 
Hence, we have performed a preliminary scanning of the parameter space, by fitting the quark 
observables in point 3 (Sect. 4.2) with ki,k^ free to float in the range -f- 3. In such a way 



we have specified three Clebsch prescriptions for k\ and k^ and they are given in eqs. (J59j, pQ| 
and (|6l|), featuring three different ansatze A, B, C, recalled here for convenience: 



ansatz A: k\ = 2 , k% = 1 

ansatz B: k\ = k^ = —3 , 

ansatz C: k\ = —3 , k% = 2 

On the other hand a priori we cannot stick to any specific value of k 2 : only the result of 
the fit will allow us to envisage the value that can give the better prediction for sin 2 29[ 2 ■ For 
that it is useful to write down explicitly^ V e2 : 

\Ve2\ 2 = cf 2 v 2 e2 + ~sf 2 v 2 eX + 2c e 12 s e 12 v e2 v el cos a (117) 

The present MSW range in eq. @) translates into \V e 2\ = (1-6 — 6) • 10~ 2 . The lower elements 
of U\ read as 

Ve2 ~ -Sl2 c 23 > v el ~ Cl2 c 13 ~ 1 > ( 118 ) 



then from eqs. ( |110[) we have |f e 2| ~ W 2 ^ ~ O.OT-^/c^ which slowly decreases with the 

asymmetry parameter 6 - for b > 8 it becomes around (5 — 6) • 10 -2 saturating by itself the 
present experimental upper limit. Therefore, the r.h.s. of eq. flll7|) tells us that the third term 
can reduce \V e 2\ if t\ 2 < 0.1 and a < 90°. In this way sin 2 26\ 2 may become smaller than 10~ 2 . 
Interestingly, it will turn out that the quark-sector fit demands the corresponding initial 12 
rotation tf 2 to lie in a similar range, if 2 < 0.1, in all the three ansatze. We also write the CKM 



elements in (87) to explicitly show the interplay between the 12 rotation tf 2 and the CP-phase 
a: 

\V US \ 2 = cf 2 vl s + ~sf 2 vl s + 2c d 12 ~s d l2 v us v cs cos a , (119) 
\Vcb\ 2 = cf 2 v 2 ch + sf 2 v 2 uh - 2cf 2 sf 2 v cb v ub cosa , 
I Kb 1 2 = cf 2 v 2 uh + sflv 2 b + 2cf 2 sf 2 v cb v ub cos a , 



(120) 



where 



Vus — C 13 S \2 ~ S Y2 i v cs — c 12 c 23 ~ ^ v cb ~ s 23 > v ub — s 13 • (1^1) 



Therefore, if the fit of the quark quantities selects a < 90° and tf 2 < 0.1, then the choice 

Je _ 2d 



k 2 = k\ is favoured as it implies t\ 2 = tf 2 . On the other hand, if it comes out that a > 90° 



29 For the sake of brevity we are assuming that the elements of the mixing matrices U e ,d are real as we have 
seen that the phases §\f are typically tiny. 
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(and irrespectively of tf 2 ) then k 2 = (i.e. tf 2 = 0) is preferable in order not to raise sin 2 26[ 2 
above 1.4 • 10~ 2 . As we shall see in the next section the three ansatze will prefer different 
values of the phase a such that in the ansatz A k% = k\ = 2 can be chosen and as a result 
sin 26\ 2 < 0.01 can be obtained. On the other hand, in both the other scenarios we have to 
select k 2 = and content ourselves with sin 2 26[ 2 > 0.01. Needless to say that this picture 
in which all the Clebsch ^1,2,3 are taken real is "minimal". For example, by allowing for 
complex k n , the lepton and quark mixing matrices would have different CP violating phases 
and therefore, for example, also in the ansatze B, C we could achieve sin 26[ 2 < 0.01 with 
non- vanishing t\ 2 . Finally, one may wonder about the alternative possibility to account for the 
large-mixing angle MSW solution of the SN anomaly requiring sin 26[ 2 ~ 0.5 - 0.98. That 
parameter range translates into \V e2 \ > 0.4 which could be achieved within our ansatze by 
large l\ 2 - l\ 2 > 0.7 - 0.8, or |a > 7 - 8. Such O(10) Clebsch coefficients seem to come less 
naturally and probably at the price of a less economical Higgs content. 



4.4 Strategy for the fit 

To study the constraints on the parameters of the model from the down-quark masses and the 
CKM mixing angles we have performed an alternative y 2- analysis. P°| 

The CKM elements \V US \, \ V c t,\,\V u b/V c b\ are assumed to be mainly affected by Gaussian 
errors and therefore the standard x 2 -function can be assumed for them: 



X 2 = E( — ) . A a = x a -x^ (122) 



where the theoretical outcome for a certain observable is denoted by x a and the corresponding 
experimental value and statistical error by x^ xp and a a , respectively (see Table Q). On the 
contrary, the uncertainties affecting the determination of the down-quark masses are due to 
the dependence on the theoretical-model used to extract the masses from the measurements. 
For this reason we think that it is not correct to assign them a Gaussian distribution and 
hence to define the canonical x 2 - function. Nevertheless, we interpret the 1-a range reported 
in the Table Q as a 'reasonable' interval accounting for all possible (different) determinations. 
Therefore, we prescribe the masses a flat distribution in that range. More precisely, we define 
the following x 2 -like function for the masses: 

2 f if |A 6 | < 0.15 

Xb = \ f\A b \-0.l5" " 



X. s 



^m^J if |Afc| > 0.15 

if |A S | < 75 

(^o^) 4 ifA s <-75 . 

^ 4 if A, > 75 



30 J 



~2 _ f if |A s/d | < 3 
Xs/d ~ 1 (\A s/d \ - 3) 4 if |A s/d | > 3 ' {US) 

30 Those observables are computed by numerically diagonalizing the Yukawa matrices ( pa) once the parameters 
A,B,D are determined from the charged-lepton masses (for given ki^,3) through eqs. (p3-Hq). The agreement 
between the numerical and the analitycal out-comings - presented in Sect. 4.1 - is at the level of one per mil. 
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Ansatz A: 


/. c = o 


II- Y ct = 


III. Complete 


IV. Spont.CP 


kx, 2 = 2,k 3 = l 




(Y t ,b,a,tj 2 ) 


(Y t ,b,£,a,tj 2 ) 


(Y t ,b,tj 2 ) 


m b [GeV] 


4.35 


4.00 


4.08 


4.04 


m s [MeV] 


219 


249 


253 


254 


ins/ " L d 


17.5 


20.8 


21.9 


22.0 


\ * US \ 


0.2138 


0.2197 


0.2193 


0.2198 


Wrh\ 


0.0507 


0.0447 


0.0437 


0.0433 


\V„h/Vrh\ 
\ v 110/ v CO \ 


0.080 


0.097 


0.094 


0.098 


★ IVmI 
A I v ta \ 


0.0148 


0.0140 


0.0133 


0.0133 


*\Vtd/V ta \ 

\ * LCI / * Lb \ 


0.30 


0.33 


0.32 


0.32 




oo 


1.9 


0.84 


0.77 


★ sin 2 20L 


0.96 


0.83 


0.87 


0.85 


★ sin 2 20[ 2 


1.2 • 10" 2 


2.0 • 10" 3 


5.2 • 10" 3 


5.7- 10" 3 


★ sin 2 26[ 3 


1.8- 10~ 2 


2.6 • 10~ 2 


2.4 • 10" 2 


2.4 • 10~ 2 


-km u /m d 


0.6 


0.4 


0.3 


0.3 


★M( /sin/? 


172.7 


172.7 


172.7 


172.7 


Y t 


0.5 


0.5 


0.5 


0.5 


b 


10.0 


12.1 


11.8 


12.0 









180° 


7T 


a 




19.7° 


43.4° 


vr/4 


id 




3.4 • 10" 2 


4.4 • 10~ 2 


4.7 • 10~ 2 


y 2 ■ & y 2 ■ 

A-min A.min 


0.06 k 50.7 


2.1 & 9.5 


0.4 & 6 


1 & 5.1 



Table 5: The analysis of the ansatz A. Taking as reference eq. ( |83| ) for Y e _d, the columns refer to: 
/. (1,3) and (3,1) entries are set to zero (C = 0), i.e. three-parameter fit Yt, b,£; II. (2,2) entry is 
set to zero (Y ct — * 0), i.e. four-parameter fit Y t ,b,a,tf 2 ', HI- complete pattern - five-parameter fit 
Y t , b, £, a, tf 2 - I n an cases, from Y t = 0.5 it comes out Y ct — 2.6 • 10~ 3 . The quantities marked by the 
symbol ★, are not included in the fits but are predictions of the fit itself in correspondence of best-fit 
parameters (also shown). The values at the minima of the (non-Gaussian) x 2 — like function Xmin ( see 
eq. ( |123| )) and of the Gaussian function, Xmin are a ^ so given. 



where x 2 , xl, X 2 s /d refer to m^, m s , m^/m^, respectively and A s , A s / d are the correspond- 
ing deviations defined as in eq. (122). We can note from ( |123| ) that mass values outside the 
l-o" range are strongly penalized by the higher-power dependence of the (arbitrary) assigned 
distribution. The sum of the above functions will be denoted as x 2 = x\ + xl + X 2 s /d- 

By minimizing the sum x 2 + X 2 ■, we rather test the consistency between the model and the 
Gaussian-distributed data within a reasonable range for the down-quark masses. In the next 
we shall report both the value of Xmin an d Xmin separately, in Tables [5L |6] and [/] containing all 
the results of the fits. We have to remark that descriptions of the quark masses outside the 
'reasonable' range may be accounted, for example, by uncertainties in a s (Mz). 
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1.0 
0.8 
0.6 



0.4 



0.2 



0.0 

-1 -0.5 p . 0.5 1 

Figure 5: The unitarity triangle as emerging from the fit of the ansatz A (case in the (p',?/) 
plane. The angles defined in eqs. ( U23 ) are to be identified as as CAB , = ABC , 7 = BCA . The 
parameters A, A are accordingly fixed as A = 0.2193 and A — 0.909. We have drawn the 1-cr contours 
of R b = 0.414 ± 0.062 (solid lines) R t = 0.913 ± 0.191 - as inferred from AM d - (dot-dashed lines) and 
Rt < 1.02 - from the lower bound on AM S (dashed- line). The band enclosed by the hyperbola (dotted) 
refers to B K = 0.8 ± 0.2. 

4.5 Testing theoretical models: fit and predictions 

We have found it to be instructive to consider first the case with (1,3), (3,1) entries set to zero 
(C = 0), to investigate the effect of a non-zero 22 entry on the pattern (|i~3|). This pattern 
(denoted by /) has three parameters, 6, £; and Y t . As a second step, we have restored the 
(1,3) (3,1) entries and set to zero the 22-entry (formally, Y c t = 0) leaving as free parameters 
b, if 2 , o", Y t (case denoted by II) . 

Finally, the complete pattern is analysed with five free parameters, b,£,,if 2 ,a,Y t (case 
denoted by III). We can first make some general considerations. We can note that the 
dependence on Y c t of the physical quantities is mainly encoded in the factors I e , Id- The 
former can be easily estimated, I e ~ 1 + 0.044 • cos £ and thus Y c t can induce quite a small 
correction on the leptonic mixing angles. Therefore, from the expressions ( |104| ), we can infer 
that |V^3| = 1 1^3 1 « s| 3 (see (pq)) depends mainly on the asymmetry parameter b, increasing 
roughly as Vb. As a result, the AN bound sin 2 26> 23 = 4|V r At3 | 2 (l — |V^ 3 | 2 ) > 0.8 requires 6 < b < 
12 as it was found in the case |^]. Analogously, in the case / it is \ V e 2\ = |f e 2| ~ s?2 c 23 
which weakly decreases with b. Then the MSW range, in the case I can be recovered for b > 7 
where sin 2 29\ 2 = 4|v e 2| 2 (l — |f e 2| 2 ) gets smaller than 1.5 • 10~ 2 . As we already know, this 
mixing can be affected by the initial 12 rotation t\ 2 present in the full model (case III). 

The discussion about 1^ affecting the quark observables is more involved as it depends, 
through the ratio Yb/(Y S — Yd) oc /c 2 , on the specific ansatz considered. However, in the case 
A a rough estimate gives Id ~ 1 + 0.12 cos £ , while larger correction can be achieved in the 
ansatze B, C, Id ~ 1 + 0.25 cos £. In the following we discuss separately our results for the 
three ansatze, reported in Tables [j| || and 0. 

Ansatz A In the case / (Table |5|- second column), the result of the fit is acceptable though, 
due to the high accuracy achieved in the determination of the CKM angles, we can observe 
some discrepancy in the Cabibbo angle |V^ S | = \v us \ w sf 2 oc Id/h and \ V c b\ ~ S23. The former 
appears to be quite small while the latter is somehow too large. These two quantities require 
in turn £ = (to maximize Id) and quite a large 23 asymmetry, b = 10 (to minimize S23) as 
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best fit points. Notice that b > 10, though preferable to further reduce \V c b\, is not tolerated 
by \V US \ which would further decrease. 

On the other hand, the quark masses fall in their 'reasonable' range. Noticeably, Yt = 0.5 
is preferred in all the fits and as a result the prediction M t j sin (3 < 172.7 GeV turns out to be 
consistent. Interestingly, the predictions (marked by *) for the leptonic mixing angle are quite 
good. We have to recall that whenever the quark mixing are 'reasonably' accommodated, the 
leptonic mixing angles fall automatically into the presently most favoured range thanks to the 
remarkable product rule (|i~2|) inherent in the theoretic structure of the model. For example 
the AN oscillatory mixing sin 2 29 l 23 is almost maximal, as the large b = 10 implies. The 13 
neutrino mixing remain well below the upper bound shown in Table 1 and sin 2 29[ 3 lies in the 
range (l-j-3) • 10~ 2 in all the ansatze, as we will see (cfr. Tables || This range in the case of 
Sm 23 close to the upper bound in (||), can be of interest for the experimental search of v e — ► v T 
oscillation in the future CERN Neutrino Factory |l0|] . 

We do not show in the tables the prediction for the neutrino mass ratios. However, as 
regards 7771/7772 the prediction is just m\/m2 = Y u /Y c (in any of the ansatze) and so it is 
typically about ~ 3 • 10~ 3 . On the other hand, the prediction for the ratio 777,2/ '^3 depends on 
the value of the asymmetry parameter b'. In the simplest case b' = b, it can be easily computed 
as 7772/7773 = bY c t and so it is in the range ~ (2 — 3) • 10~ 2 in agreement with the present 
experimental hint. This out-coming will be similar in all the ansatze and it is a byproduct of 
the link between the neutrino and the up-quark masses and the large 23-asymmetry b demanded 
by the quark phenomenology. 

The other predictions, Vtd, Vtd/Vt s are marginally compatible with the experimental values. 
We have also reported the expected mass ratio m u /md extracted from the 'ellipse' parameter 
Q (116) which appears to be consistent with the determination of m^/m^ obtained from the 
fit. 

Clearly there is no CP violation in the CKM matrix as the infinite value of Bk does reflect. 
The situation gets improved in the case II. The presence of the (1,3), (3,1) entries introduces 
the initial 12 rotation tf 2 and the phase a which strongly modifies the CKM elements (|S7|). 
From the expressions of the Cabibbo angle in ( |119| ) and from the fact that the element v us 
itself is around 0.2, we can deduce that Of 2 is to be a small rotation, t\ 2 ^0.1, to prevent 
too large a correction from the element v cs ~ 1. So \V US \ is increased up to its experimental 
value by tf 2 ~ 0-03 and a ~ 20°. At the same time \V c f,\ is reduced, thanks to the larger 
value of the asymmetry parameter, b ~ 12 which is now not prevented by V us . We can notice 
that also (V^/V^l is successfully reproduced. The effect of b larger, i.e. of the decreasing 
of the ratio Y^/Y T has in turn induced a bigger m s and a smaller rrif,, both slightly outside 
their 'reasonable' range. The amount of CP violation from the CKM matrix, though increased 
thanks to the non- vanishing phase a, is still not enough, Bk ~ 2. The predicted values of 
the leptonic mixings are good. Accordingly to the approach elucidated in the Sect. 4.4, we 
have fixed the initial 12 lepton rotation as t\ 2 = tf 2 , i.e. = k\. Then, as expected, sin 2 29[ 2 
is strongly reduced below 10 -2 with respect to the case /. Notice, that the amount of CP 
violation in both the quark and lepton mixing matrix is the same as it is controlled by the 
same phase a. 

Finally, by comparing the results obtained in the case III - where the 22-entry is restored 
- we conclude that the presence of the 22-entry does not play a significative role - the quality 
of the fit is quite stable in the two cases. However, the predictions are better. Notice that 
the right amount of CP violation is achieved - Bk ~ 0.7 - as the phase a is now bigger. For 
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Ansatz B 


J. C = 


II. Y ct = 


Complete 


IV. Spont.CP 


&i,3 = —3, k 2 = 




(Y t ,b,a,tj 2 ) 


(Y t ,b,£,a,tf 2 ) 


(Y t ,b,tj 2 ) 


m b [GeV] 


4.65 


4.82 


4.24 


4.37 


m s [MeV] 


76 


83 


87 


83 




20.3 


27.3 


24.4 


23.0 


| v US | 


0.2127 


0.2194 


0.2196 


0.2194 


\Vrh\ 
\ r C0\ 


0.0430 


0.0575 


0.0403 


0.0409 


\ v 110/ v CO \ 


0.014 


0.102 


0.094 


0.099 


* \ v ta\ 


0.0097 


0.011 


0.0086 


0.0089 


*\Vtd/V t8 \ 

\ r LLL / " Lb \ 


0.23 


0.20 


0.22 


0.22 




oo 


0.22 


0.69 


0.62 


★ sin 2 28L 


0.97 


0.64 


1 


1 


★ sin 2 2d[ 2 


1.5 • 10~ 2 


1.7- 10" 2 


1.3 • 10~ 2 


1.4 • 10" 2 


★ sin 2 29[ 3 


1.0 • 10~ 2 


4.5 • 10~ 3 


1.5 • 10~ 2 


1.4 • 10" 2 


-km u /m d 


0.4 


— 


— 





*M t f sin (3 


172.7 


195.6 


172.7 


172.7 


Y t 


0.5 


0.97 


0.5 


0.5 


b 


6.8 


3.5 


8.8 


8.2 









9.7° 





a 




95.7° 


93.5° 


vr/2 


z 12 




0.1 


9.2 • 10" 2 


9.7 • 10" 2 


y 2 ■ & y 2 ■ 

A.min A.min 


11.4 & 45.2 


150.5 & 112 


0.02 & 0.24 


& 0.82 



Table 6: The analysis of the ansatz B. In the fit II from Y t = 0.97 it follows Y ct = 1.9 • 1CT 3 . See also 
the caption of Table B. 



the same reason, also sin 2 26 l 12 lies exactly in the range required by the MSW solution. 

The fact that the best-fit points of the phases, £ = ir, a ~ vr/4, are just integers (or half- 
integers) of 7r may suggest that the CP violating phase is originated by a spontaneous-breaking 
mechanism. For this reason in the column IV we have performed again the fit fixing the phases 
£ and a at the nearest "7r-fold" values of the kind qn (q = 1, 1/2, 1/4, • • •) as respect to those 
given by the fit III. In this case we do not observe any variation. Therefore we can conclude 
that the ansatz A provides quite a good fit of the quark observables for large 23 asymmetry 
b ~ 12 and tf 2 ~ 0.05 and for £ = 180° , a = 45°. The corresponding amount of CP violation 
appears to be in the right range as indicated by Bk ~ 0.8. We shall come back to the CP 
issue below. 

Ansatz B We have first to remind that this Clebsch-pattern in its Fritzsch version (vanishing 
22 entry, C = and 6=1) would give rris/m^ ~ 25 and hence too a small | V^ s | ~ 0.20. Let 
us consider the effect of the non-vanishing 22 entry. From eq. ( |102[ ) we see that a proper 
trend can be achieved for £ ~ and b > 6 so that to maximise Id- On the other hand, the 
moderate asymmetry b ~ 6 — 7 required by the fit implies quite a large nib, small m s and 
consequently also a very tiny \Vub/Vcb |. The fit is indeed rather poor (see Table |(| column 
/). However, setting to zero the 22-entry and taking C ^ spoils completely the fit (Table 



40 



1 


1 1 


1 


1 


i j, ■ i | i _j, i 


i . 










— 




- 
- 


(B) 




/ 

/ 


/ • ' ^ 


; - 
.'— * ~ 


- 
- 






/ 

i 


/ • ' ' 

/ 

/.-'' / 


- 














- 






T / 

/i ' 7 












/-■■ I 








i i 


i llil 


■ ' 

' 1 1 
1 , ll/c 


i l l 

i i i i u i i i 




1 




-0.5 


p- 


0.5 


1 



Figure 6: As in Fig. (§) for the ansatzB (case III). Here A = 0.2196 and A = 0.836, R b = 0.413±0.062, 
Rt = 0.989 ± 0.207 (from AM d ) and R t < 1.02 (from AM,). 

||, column II). Among the CKM elements only \V US \ and the ratio \V u b/V c t>\ get improved 
thanks to the tf 2 ~ 0.1 rotation. The same rotation is instead less important for \V c b\ which 
remains too large due to the small asymmetry, b ~ 3 — 4. For the same reason also m s /md 
is quite large. In particular as m s /rrid is larger than the parameter Q, the light-quark mass 
ratio m u /md cannot consistently be evaluated. As a whole this case is definitely disfavoured. 
Remarkably, the leptonic mixing angles as well as Bk are not well predicted, too. 

Nevertheless, the interplay of the 22-entry with the (1,3), (3,1) entries can offer a satis- 
factory description as the results in the case III show. All the quark quantities are perfectly 
fitted (notice that both Xmin an d Xmin are < -Q- The fit requires b ~ 9, tf 2 ~ 0.1, small £ and 
large CP phase, a ~ 94°. All the predicted quantities are within their experimental ranges. 
The 23 leptonic mixing is predicted to be maximal as required by the AN anomaly. The large 
phase a enforces if 2 = or k 2 = and the SN oscillation mixing comes out to be close to 
the upper limit, sin 2 29[ 2 = 1.3 • 10~ 2 . Indeed, were t\ 2 = tf 2 (k 2 = k\) and a ~ 90° - as 
preferred by the fit - the 12 leptonic mixing |V^| would get further increased (see eq. ( |11T[) ), 
contrary to what happens in the ansatz A. As a further consequence of t\ 2 = the CP phase 
is vanishing in the lepton mixing matrix except for the small contribution driven by the phase 
£. Finally in the last column, we have considered the case with £ = and a = tt/2 fixed. The 
quality of the fit remains very good. 



Ansatz C The Fritzsch case for this pattern shows an opposite behaviour as compared to 
the previous ansatz. Indeed, it is m s /md ~ 12 implying too a big | | ~ 0.29. Therefore, 
in the variant /, this can be cured with £ > 90° and a large b (see Table 0, column I). In 
this case only iVub/V^I is strongly incompatible with its experimental range. Notice that, for 
the same asymmetry b = 10, m;, is smaller as respect to the value obtained in the ansatz A 
(I), since as k\ is larger as Z becomes smaller. All the predictions are reasonable (except for 
Bk). In the second case (II) all the CKM mixings are well reproduced thanks to the initial 
12 rotation, but Y c t = makes m s /m^ below the 'reasonable' range. Finally, the complete 
model provides quite an excellent fit (III)- Interestingly, also this ansatzpoints to maximal CP 
violation - a ~ 90°. In correspondence of the best-fit parameters all the predicted quantities 
are in agreement with the present experimental status. We should remark that like in the 
previous ansatz, the presence of the 22-entry is crucial. Moreover, also in this case the choice 
k 2 = does not induce CP violation in the lepton mixing matrix. In the last column, we 
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Ansatz C: 


L C = 


II. Y ct = 


III. Complete 


IV. Spont.CP 


h = -3, k 3 = 2, k 2 = 




(Y t ,b,a,tj 2 ) 


(Y t ,b,t,a,i( 2 ) 


(Y t ,b,tj 2 ) 


m b [GeV] 


4.04 


4.02 


4.17 


4.09 


m s [MeV] 


115 


111 


114 


120 


nog/ i * i>a 


18.2 


16.7 


18.9 


20.4 


\ Y us \ 


0.2197 


0.2196 


0.2195 


0.2205 


\Vrh\ 
1 CO | 


0.0378 


0.0381 


0.0395 


0.0381 


| v 111) J v CO \ 


0.041 


0.093 


0.093 


0.082 


A I v ca | 


0.0098 


0.010 


0.0098 


0.0093 


★ IKd/VLl 

| * tti / * Lj i 


0.27 


0.28 


0.26 


0.25 




9.3 


0.74 


0.65 


0.83 


*sin 2 2^o 


0.97 


0.97 


1 


0.99 


★ sin 2 2^ 2 


1.2 • 10" 2 


1.2 • 10" 2 


1.3 • 10~ 2 


1.3 • 10~ 2 


★ sin 2 20^ 


1.8 • 10~ 2 


1.8 • 10~ 2 


1.4 • 10~ 2 


1.6 • 10~ 2 


-km u /m d 


0.6 


0.7 


0.5 


0.4 


*M t / sin /3 


172.7 


172.7 


187.0 


181.2 




0.5 


0.5 


0.7 


0.6 




10.0 


10.1 


8.7 


9.4 




107° 




135° 


7T 






106° 


93.3° 


vr/2 






7.3 • 10" 2 


8.6 • 10" 2 


7.3 • 10" 2 


Y 2 ■ & y 2 ■ 

A mm A, mm 


0.5 & 15 


4.0 & 0.7 


& < 0.01 


0.02 & 1.5 



Table 7: The analysis of the ansatz C. In the fit III (Y t = 0.7) and IV (Y t = 0.6) it is Y ct = 2.2 ■ 10~ 3 
and 2.4 • 10 -3 , respectively. See also the caption of Table (||). 



have performed a three-parameter fit, setting a = 90° and £ = 180°. One could conclude that, 
by endowing the same pattern with spontaneous breaking of CP-symmetry, the description 
remains very good. 

Finally, we can try to give an overlook at all the three ansatze analysed. The case A requires 
the CP-phase a < 90°, whereas both the ansatze B and C prefer maximal CP-violation, 
a > 90°. This should imply in general a different relations between the corresponding CKM 
elements. For example the three ansatze will show up a different shape of the "unitarity" 
triangle characterized by the following angles 

Therefore in correspondence of the results of the fits obtained in each ansatz for the com- 
plete case III, we have plotted the unitarity triangle in the p',rj' plane (Fig. 6, 7, 8). We 
recall that p,r/,A,\ are the Wolfenstein parameters and p 1 = p(l — A 2 /2), rf = 7/(1 — A 2 /2). 
More precisely, the side CB corresponding to V c dV* b has been chosen real and rescaled to unit 
length. Hence the length CA and CB in the rescaled triangle, usually denoted by R b and Rt, 
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0.8 



0.6 — 



0.4 



0.2 



0.0 




Figure 7: As in Fig. (§) for the ansatzC (case III). Here A = 0.2195 and A = 0.820, R b = 0.413±0.062, 
R t = 1.01 ± 0.21 (from AM d ) and R t < 1.02 (from AM,). 

respectively, are: 



Rb = 
Rt = 



v ud v 



ub 



VcdVZ 

VtdVl 



pi2 _|_ 



v 2 j A 



(l-p') 2 +ri 



12 



1 



Vcb 



Once the parameters A, A have been fixed according to the results of the III fit, A = | V u 



(125) 

(126) 

,A = 
(solid 



l^usl/A 2 , we have depicted the 1-a circles of Rb from the measured value of 
lines), and the circles of Rt from the the mass difference AM^s decribing the strength of 
the B® s — B® s mixings (dot-dashed and dashed lines, respectively). The values used are 
AM& = 0.471±0.016ps _1 and AAf s > 12.4ps _1 [41]. In the same plane the hyperbola represent 
the iso-contours for Bk (dotted lines). This overall picture of the CKM elements shows that 
the ansatzA (see Fig. |5|) is characterized by a large angle 7, 7 > 100°, and by a small angle 
j3, (3^8°. They both appear to be in disagreement with the result reported by the CDF 
collaboration sin 2(3 = 0. 79^44 and with the indirect bounds on 7 which disfavour 7 > 100° 
Moreover, the vertex A lies above the region delimited by the circles Rt due to the large 
value of \Vtd\ ~ 0.013 predicted. On the contrary, the ansatzB provides a consistent scenario 
where the observed amount of CP-violation can be accounted by the CKM picture (Fig. |6|). 
The same holds for the ansatzC (Fig. ^). Ultimately, a better discrimination of the specific 
ansatz will be possible when the program of the measurements of the present and future B 
factories is completed. 



5 Conclusions 

We have considered a particular set of textures for the Yukawa matrices of quarks and leptons 
which share some features with both the Stech ansatz proposed long time ago Jl6|, [l7]] and the 
Fritzsch-like textures suggested in [13]. Our main motivation was to show that they are not 
only successful from the phenomenological point of view but also grounded and compelling 
on the theoretical side. We have discussed how these textures could emerge in the context of 
grand unified theories, in terms of the prototype SU(5) model complemented by the horizontal 
symmetry SU(3) H (Sect. 2). 
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The SU(3)h group may seem too 'large' as compared to U(2)h- Indeed the latter, providing 
ab initio the 2 + 1 representation structure for the fermion fields, directly singles out the heavy 
top quark from the lighter fermions of the first and second generation. Namely, the models 
based on U(2)h J4*| invoke the familiar paradigm according to which the third generation has 
a priori order 1 (tree-level) Yukawa couplings while the small Yukawa constants of the lighter 
generations emerge from higher-order terms containing the U(2)h symmetry breaking Higgses, 
with VEVs smaller than the cutoff scale. 

However, there are many good points in favour of SU(3)h- First, it accounts by itself for 
three fermion families, and thus can be more predictive than U(2)h- Second, the spontaneous 
breaking features of SU(3)h may turn the Yukawa constants of the low energy theory (MSSM) 
into dynamical degrees of freedom and fix the inter-family hierarchy in a pretty natural manner. 
Namely, the third generation becomes heavy (Yt ~ 1), while the second and first ones become 
lighter by successively increasing powers of small parameters. 

In the general discussion of theoretical issues, we have put forward the following argument: 
just like the SU (5) theory contains the adjoint Higgs, 24-plet, which, by breaking down SU(5) 
to SU(3) x SU(2) x U(l), does hide the existence of the large GUT symmetry, also SU(3)h 
may be first broken down to U(2)h by its adjoint Higgs, an octet, which makes the existence 
of a larger symmetry among all three families less visible. 

In particular, the horizontal octet which breaks SU(3)h down to SU(2)h, not only diver- 
sifies the third family from the lighter ones, but also induces a non-trivial Clebsch structure 
between the 1-2 and 3 generations in the Yukawa matrices of quarks and leptons. The effect of 
the latter is twofold. First, it induces the complementary "large - small" (see-saw like) pattern 
between the neutrino and quark mixing angles, which is indeed exhibited by the observed small 
value of V c b and the nearly maximal j/„ — v T mixing. Second, it may link the neutrino mass 
hierarchy to that of the up-quarks in realistic way. 

On the other hand, the SU(5) adjoint has proved necessary to break the 'standard' down- 
quark and lepton degeneracy, providing different Clebsch coefficients in the Yukawa matrices. 
Though its role is not univocal, we have featured three Yukawa patterns characterized by what 
we consider as the most 'natural' Clebsch factors. In Sect. 3 we have elaborated a consistent 
50(10) x SU {3)h model which, supplemented by the symmetry Q = U {I) a xZqxTZ, represents 
an existence proof of the Yukawa texture suggested. In particular, the specific model outlined 
provides an understanding of the Clebsch pattern of the ansatz B. Remarkably, the same 
symmetry content of Q also motivates the missing VEV mechanism necessary to achieve the 
doublet-triplet splitting. 

A careful fit of the quark observables demonstrates that all the ansatze A, B and C 
reproduce quite well the down-quark masses and CKM mixing angles (Sect. 4). What appears 
very interesting is that for the best-fit ranges of the parameters the predictions for the neutrino 
mixing angles are in good agreement with the present experimental hints. About the lepton 
mixing angles, we have to stress that while the outcome for the 23 neutrino mixing is a 'genuine' 
prediction, in strict connection only with the down-quark sector through the product rule (12), 
the 12 mixing angle becomes effectively a free parameter fixed by the rotation t\ 2 - Should 
the small angle MSW solution be confirmed in the future, then all the ansatze can account 
for that with different values of t\ 2 - On the contrary, should the 12 neutrino mixing angle be 
much larger, as indicated by the large-mixing angle regime of the MSW solution, then only 
the lepton sector would need some refinement. 

Interestingly, also the prediction for CP violation in the quark sector may lie in the observed 
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range. In particular, all the ansatze generate a CKM matrix with the correct amount of CP 
violation. Though the experimental test of CP violation in the SM is not yet accomplished, this 
result may imply that superpartners contributions to ek must be adequately suppressed pql 
and it would be interesting to study in the theoretical models presented here the implications 
of the horizontal symmetry on the SUSY spectrum. 
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Appendix A 

Consider the following superpotential terms including antisextet superfield S = S 13 and its 
sextet partner S = Sif. 

W S = -^SS + S 3 + S 3 , (127) 

where, to be explicit, S 3 = \eijk£imnS l1 S^ m S kn = det S (similarly for S 3 ), and order 1 coupling 
constants are absorbed. Observe that this superpotential is manifestly invariant under Z3 
symmetry: S — > ex.p(i^-)S and S — > exp(— i^-)S. 

Without loss of generality, the VEV of S can be chosen in the diagonal form, (S) = 
Diag(5i, 52, 53). Then the condition of vanishing F-terms F$,F§ = implies that (S) is also 
diagonal, (S) = Diag(5i, 52, 5?), and 

5152 = //53, 5i52 = /i53, 

5153 = /i5 2 , 5i53 = ^5 2 , 

5 2 5 3 = //5 1 , 5 2 5 3 =fiSi. (128) 

In the exact super symmetric limit the VEV pattern of S and S is not fixed unambiguously and 
there remain flat directions which represent a two-parameter vacuum valley. In other words, 
the six equations fll2S| ) reduce to four conditions: 

5i5i = 5 2 5 2 = 5 3 5 3 = /u 2 , S 1 S 2 S 3 = fi 3 ; S^S^fi 3 , (129) 



while the others are trivially fulfilled (the last in eq. (|129j ) is not independent, it just follows 
from the others). Thus, in principle the eigenvalues 5i j2j 3 can be different from each other, 
say 53 > 5 2 > S\. Then eqs. Ql28p imply that 5i i2j 3 should have an inverse hierarchy, 
53 < 5 2 < Si. More precisely, we have Si : S 2 ■ S3 = S^ 1 : S^ 1 ■ S^ 1 . 

The flat directions of the VEVs can be fixed from the soft supersymmetry breaking terms. 
The relevant ones are D-termsP 1 ] 



£ = - f d*9zz aTtS^S + -^(TrS^) 2 + -^Tr(S^SS^S) + ... 
J [_ M z M z 



+ [S -> S] (130) 



31 The F-terms~ f d 2 8zW are not relevant for the VEV orientation as far as they just repeat the holomorphic 
invariants like SS and det S whose values are already fixed by the conditions ([129) . 
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am 2 (|5 1 | 2 + |5 2 | 2 + |53| 2 )+/?^(|5 1 |V|5 2 | 2 + |53|Y+7^(|5i| 4 +|5 2 | 4 + |5 3 | 4 ) + [|5 fc |-|5 fc 



having a similar form for S and S. Here z = m6 2 and z = mO 2 are supersymmetry breaking 
spurions, with m ~ 1 TeV. The cutoff scale M is taken as the flavour scale, i.e. the same as 
the one in superpotential terms (|24|), and we assume that M 3> fi. Therefore, for Si 23 these 
terms translate explicitly into the following scalar potential: 

A ; ( i^+i^+is3iv+7^ 

(131) 

The stability of the potential implies that (5 > and 7 > —/3, whereas a can be positive or 
negative. In the former case the minimization of the potential ( |131| ), under the conditions 
( |129| ), would imply that Si = 5 2 = £3 = fi, i.e. no hierarchy between the fermion families. In 
the latter case, however, the largest eigenvalue of S and S, respectively £3 and Si, grow up 
above the typical VEV size \x and reach values of the order of the cutoff scale M: 

s ^ = (,W^f M < 132 > 

that is S3, Si ~ M. Then it follows from (|129|) that 



S 2 ,S 2 = n = eS 3 , Si,S 3 = — = e 2 S 3 , (133) 

<->3 

where e = ///5 3 ~ ji/M. 

Let us assume now that the theory contains three triplets superfields A n = A n i, and their 
partners A n = A l n , n = 1,2,3. One can incorporate them by the following terms in the 
superpotential: 

W A = -J2 HnA n A n + ^i^2^3 + ^1^3 (134) 

n 

where order 1 coupling constants are understood and A1A2A3 = E^AuAzjAsk. For simplicity, 
we shall take all masses \i n equal, ^1,2,3 = yu' < M. In addition, by assuming that the triplets 
have Z 3 charges different from that of S, we do not include terms like SA 2 = S ij ' A^Aj. In 
this way, in the exact supersymmetric limit the ground state has a continuous degeneracy (flat 
direction) related to unitary transformations A n — > UA n with U C SU(3)a- In other terms, 
the superpotential W = Ws + Wa has an accidental global symmetry SU(3)s x SU(3)a, with 
the two SU (3) factors independently transforming the two sets of horizontal superfields, S and 
A. 

Similarly to the case of the Higgs fields S,S, now the conditions Fa, Fa = can only fix 
the values of holomorphic invariants A n A n and A1A2A3. Namely, by unitary transformation 
A n — > UA n (U C SU(3)a), one can choose a basis where the VEV of Ai points towards the 
first component. Then we see that in this basis the fields A 2 and A3 should have the VEVs 
towards the second and third components. In other words, in this basis the triplets have the 
VEVs (A n i) = 5 l n Ai, satisfying the following equations: 

A1A1 = A 2 A 2 = A3A3 = ii 2 , ^i^ 2 ^3 = /"' 3 , A1MA3 = /d' 3 . (135) 

The three VEVs can be different and one can take, say Ai > A2 > A3. The hierarchy between 
the latter can be fixed by soft D-terms 



C = - d 4 6zz 



a ' ^ ^-n^n + ~TP} y^X^ln^rn)(Al l A n ) + ^ {A\ n A n ){A\ l A m ) 

n lV1 ran 1V1 m ^n 



~t~[/4n > A n ] 

(136) 
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with a' < 0. Then, similarly to what happened to the sextets, the largest VEV rises up to 
the cutoff scale M: Ai,A 3 ~ M, while A 2 ,A 2 ~ e'M and A 3 ,Ai ~ e /2 M, with e' ~ fj//M ~ 

As for the relative orientation of the sextet and triplet VEVs, these will also be fixed from 
the soft D-terms: 

C = - [ d A 0zz-^ £(4SSU n ) + [ S _ S, A n -» I n ] (137) 

Indeed, for positive <5 the favoured orientation between the <S and A bases, which minimizes 
the energy of the ground state, corresponds to U = 1. 

The above can be interpreted in the following manner. In the view of operators like (p4[), 
the MSSM Yukawa constants become dynamical degrees of freedom. In particular, the operator 
O u in <M) implies 



( Y u i 

Y c ] - K - T7 L = i 7 & |. (J.3S) 

y o o 



M 




In the exact super symmetric limit the values of the constants Y UCt t are not fixed - they have 
flat directions and eq. ( |129| ) translates into the constraint Y u Y c Yt = e 3 . However, we have 
shown that soft supersymmetry breaking terms could naturally split the Yukawa constants 
so that Yt ~ 1, while Y c ~ e and Y u ~ e 2 , which perfectly reflects the observed pattern if 
e ~ 5 • 10~ 3 . (For discussions on dynamical Yj in supersymmetry see also ref. [43].) 

As for the triplet VEVs, in view of the operator O (^4|) they provide off-diagonal Yukawa 
entries 

/ A 3 A 2 \ Mi n i / A 3 A 2 \ 

A=\-A 3 A X \-Y, { -W = m\-- ^3 A 1 \. (139) 
\ -A 2 -At J n y -A 2 -A! J 

with the constraint A\A 2 A 3 = e' 3 in the supersymmetric limit, while the soft breaking terms 
can dynamically fix their values as A\ ~ 1, A 2 ~ e' and A3 ~ e' 2 . In the context of our work 
it is clear that such a hierarchy of the off-diagonal entries in Y^e is also well suited for the 
observed fermion mass and mixing pattern if e' ~ 0.1. 

One has to remark that the above considerations are literally valid if the horizontal sym- 
metry is global. For the case of local SU(3)h , apart from F-flatness conditions, there are eight 
additional conditions: the gauge D terms D a = J2n X^T a X, should vanish on the vacuum 
configuration (T a are SU(3)h generators, a = 1,..8). Although this does not occur for the 
obtained VEV pattern of S and A n , one can easily imagine the theory to contain additional 
"spectator" superfields in different representations of SU(3)h and their VEVs are oriented so 
that to cancel the contributions of (S) and (^4) in gauge D-terms of SU(3)h- In order not to 
affect the obtained solutions for S and A n , these extra superfields should not couple to the 
latter in the superpotential. 



Appendix B 

Let us assume now that the considered theory has an exact CP invariance. In other 
words, we can choose a superfield basis where all coupling constants in the theory are real, 
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including the coupling constants with the matter superfields, i.e. those in the terms (^4|). 
Applied to the superpotential (|127| ) this means that the dimensional parameter /i as well as 
the understood constants at 5 3 and S 3 are real. Without losse of generality, by means of a 
SU(3)h transformation, one can chose a basis where £3 and 1S2 are real and positive. Then, 
since eq. ( |129| ) tells us that the VEV products S1S2S3 = fi 3 is real, also the component 5i 
should be real. The same consideration is true for the VEVs of triplets A n . Therefore, no 
CP-violation can occur in this case. 

Let us now consider an alternative superpotential for triplets ( |134| ) not containing the mass 
terms but invoking instead an auxiliary singlet I: 

W A = K n IA n A n - A X A 2 A Z - A X A 2 A Z + 3A 2 / + J 3 , (140) 

where A is a dimensional parameter, K n are coupling constants and order one constants are 
also understood in other terms. Applied to this superpotential, CP invariance means that K n 
as well as A 2 are real. For simplicity, let us assume again that all n n are equal, K, n = k, and 
k < 1. 

It is clear, that the i^-term conditions Fa,F^ = are the same as in ( |135| ) apart from the 
fact that the mass scale // should be substituted by ttZ, where X = (I). The latter is then 
fixed by the condition Fj = given as 

1 3 

A 2 + J 2 + -« A kA k = A 2 + (1 + k 3 )X 2 = (141) 
3 fc=i 

Hence, X 2 ~ —A 2 , and if A 2 is positive, then X should be imaginary and so we obtain spon- 
taneous CP violation. In other words, imaginary X means that the CP-odd component of the 
superfield I acquires a VEV. 

As a result, the 'induced' mass // = kX and so the VEV product A\A 2 A 3 = // 3 are also 
imaginary. Therefore, at least one of the VEVs -4-1,2,3 should be complex and thus would 
induce the corresponding phase in the fermion Yukawa matrices. 

To fix the point, we have to find the relative phase orientation between the VEVs of S and 
A. These should be determined by soft terms like 

C = - J d 4 ezzj^(S^SU n A n + h.c), (142) 

etc., with £ being a real constant. In particular, the above term induces the following couplings 
in the Higgs potential: 

2 

C^(S 3 S 2 Al + SsSiAl + S 2 S 1 A 2 3 ) + h.c. (143) 

We see that for £ > the minimization of the ground state energy requires that in the basis 
where Sk are real and positive, A\ 2 3 < 0, i.e. all triplets have imaginary VEVs. For ( < 
we see instead that Ai i2 prefer to be real. From the couplings (143) also the smallest VEV ^3 



would prefer to be real; however it is forced to have a complex phase ir/2 by the constraint that 
the product ^4i^42-^3 is imaginary. Both these cases can be of interest for the CP-breaking 
phases of the fermion Yukawa textures discussed in this work. 
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